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DUST EMISSION AND STAR FORMATION IN STEPHAN'S QUINTET 
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ABSTRACT 

We analyse a comprehensive set of MIR/FIR observations of Stephan's Quintet (SQ), taken with the 
Spitzer Space Observatory. Our study reveals the presence of a luminous {Lm « 4.6 x 10'''^ erg/s) 
and extended component of infrared dust emission, not connected with the main bodies of the galaxies, 
but roughly coincident with the X-ray halo of the group. We fitted the inferred dust emission spectral 
energy distribution of this extended source and the other main infrared emission components of SQ, 
including the intergalactic shock, to elucidate the mechanisms powering the dust and PAH emission, 
Q I taking into account collisional heating by the plasma and heating through UV and optical photons. 

Combining the inferred direct and dust-processed UV emission to estimate the star formation rate 
(SFR) for each source we obtain a total SFR for SQ of 7.5 M0/yr, similar to that expected for non- 
I interacting galaxies with stellar mass comparable to the SQ galaxies. Although star formation in SQ is 

mainly occurring at, or external to the periphery of the galaxies, the relation of SFR per unit physical 
area to gas column density for the brightest sources is similar to that seen for star-formation regions in 
galactic disks. We also show that available sources of dust in the group halo can provide enough dust to 
produce up to Lju « 10^^ erg/s powered by collisional heating. Though a minority of the total infrared 
(-H ' emission (which we infer to trace distributed star-formation) , this is several times higher than the X-ray 

Q^. luminosity of the halo, so could indicate an important cooling mechanism for the hot IGM and account 

for the overall correspondence between FIR and X-ray emission. We investigate two potential modes of 
star-formation in SQ consistent with the data, fuelled either by gas from a virialised hot intergalactic 
medium continuously accreting onto the group, whose cooling is enhanced by grains injected from an 
^ ' in-situ population of intermediate mass stars, or by interstellar gas stripped from the galaxies. The 

former mode offers a natural explanation for the observed baryon deficiency in the IGM of SQ as well 
as for the steep Lx-Tx relation of groups such as SQ with lower velocity dispersions. 

^ I Subject headings: galaxies: groups: individual (Stephan's Quintet) - galaxies: evolution - interactions 

' - ISM - intergalactic medium 

(N ■ 

. 1. INTRODUCTION the intergalactic medium (IGM) is being converted into 

T,, . , . . . , r stars. Whereas at early epochs gas fuelling of galaxies 

. . Physical processes occurrmg m the environments oi . , , , , , i • i i ■ • i 

O r 1 • , r J i 1 1 • J J- • IS thought to proceed via cold gas accretion m low mass 

' groups ot galaxies play a fundamental role m determm- t^h^ttt ,i • i • • i- , i , i • i -i •■ i i ,i 

^ ■ inp- the star formation histnrv of the Universe GaJaxv DMHs, this mechanism IS predicted to be inhibited by the 

' ' , my iiie biai loimaiioii iiioLoiy oi iiie uiiiveibe. vjaiaxy i • i • • i , , r ,i x/-^t\ r • ^ • i t-vtv ttx 

O. , 1 . 1 J. , 1 1 , higher virial temperature ot the IGM m high mass DMHs 

■ groups are associated with intermediate mass dark mat- , ,. , t~, — ; — nr^^^^\ ,i 

^— I . Z 1 1 ^r^n^tJ^ u- u • i 1 -J-- • hostmg galaxy groups (Dekel et al. 2006 . I^urthermore, 

' ter haloes DMH which occupy a pivotal position m the , r ,i i • i , T"^ i i i 

r i- r i i it because ot the high densities, galaxy-galaxy and galaxy- 

tormation ot structures, acting as a centre ot aggrega- t^h^ • , i i, i rr- • ■ • 

,. r, TMVTTT J -Li ■ • J- J 1 ■ IGM interactions should become ertective m removing m- 

, J tion ot lower mass DMHs and their associated galaxies , , „ r ^ ■ ■ i 

u-i 1^ • ii, 1^ -ij- ui 1 r j-i i -1 terstellar gas from galaxies m groups, leading ultimately to 

rS while being the building blocks tor the mos t massive cms- , r . % ,• ? , i i • ,• i 

Lj', £i-ii,ir \. Ha — '• i — 7 — n a quenching ot star formation m the already existing galax- 

!3 . ters ot galaxies that form at later epochs (ISprmgel et al.l ■ , ■ , i r n ■ , ,i ^ ,i ,i i ^ 

^ . lonnch T ii, i. ■ u J- cn ^ fTi Ti les which have lallen into the groups. On the other hand, 

- - ' 120051 ). In the present universe about 50 percent oi all stel- ,, . , ,. , , , , . , . , , 

, . , • J -ii • r i J- 1 1-1, these same interactions lead to a chemical enrichment as 

lar mass is contained withm groups ot total mass higher t^ht i ■ i -,1 ,i , • ^ i ■ tm^t 

, 71^/11^^ ^loc rcT~^ — nionncK the IGM bccomes mixed with the stripped galaxian ISMs 
than loamM/Mr^. > 12.5 flEke et al.ll2005l ). The group en- ... , ^. , ..^^ ^ ^, 

, rr i 7i J i-i. • 1 • 1 1- • • which, potentially, could enhance the cooling ot the IGM 
vironment anects the modality m which baryomc gas m j r- jj o 

^ Max Planck Institute fiir Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, Germany; giovanni.natale@mpi-hd.mpg.de, richard.tuffs@mpi- 
hd. mpg.de 

^ Infrared Processing and Analysis Center, California Institute of Technology 100-22, Pasadena, CA 91125, USA 
^ University of Central Lancashire, Preston, PRl 2HE, UK 

* Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 Saint George Street, Toronto, ON, M5S 3H8, Canada. 
^ Department de Fi'sica Teorica y del Cosmos, Universidad de Granada, Granada, Spain 

^ NASA Herschel Science Center, IPAC, California Institute of Technology, Pasadena, CA 91125, USA 

^ Research School of Astronomy & Astrophysics, The Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia 

* Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, Dapnia/Service dAstrophysique, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, 
France 

^ Purple Mountain Observatory, Chinese Academy of Sciences, 2 West Beijing Road, Nanjing 210008, China 
^° Spitzer Science Center, IPAC, CaUfornia Institute of Technology, Pasadena, CA 91125, USA 

Instituto de Astrofi'sica de Andalucia, CSIC, Apdo. 3004, 18080, Granada, Spain 
^■^ Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 

1 



2 



Nat ale et al. 



and its ability to accrete onto existing galaxies and to form 
new star forming systems. The relative importance of all 
these processes in determining the star formation activity 
which is observed to be taking place in groups is an open 
question. 

The Stephan's Quintet compact group of galaxies (SQ) 
[3 presents a natural laboratory with which all these phe- 
nomena affecting star form ation in g r oups can be stud- 
ied. As shown bv .Trinchieri et al.l (|2005D (T05) and 
lO'Sullivan et all (l2009|) (0S9), this group presents a dif- 
fuse halo of X-ray emission extending in radius at least as 
far as 40 kpc with the bulk of the gas radiating at temper- 
atures « 6 X lO^K. Assuming hydrostatic equilibrium, this 
indicates a dark matter halo mass of w IO^^Mq, interme- 
diate between galaxies and clusters. The metal abundance 
of this hot gas is rather poorly constrained, consistent with 
a primordial and/or galaxian origin. 

One galaxy, NGC 7318b, is apparently unbound, enter- 
ing the group at a high relative velocity of « 1000 km/s 
and colliding with the group IGM, as evidenced by 
a ^ 32 kpc no rth-south long ridge prominent in ra- 
dio continuum ()Xu et al.l 12003.), optical line emission 
(|Xu et al.l IT999I) X-rav (|Trinchieri et al.l l2005l) and re- 
cently also powerful m i d-infrared rotational hydrogen lines 
(jAppleton et al.l l2006l : ICluver et all I2OIOO . The galax- 
ies of the group present extended tidal tails that have 
been used to constrain their recent i nteraction history 
(|Moles et al.lll998l ISulentic et al.ll2001t) . Potential conse- 
quences of these interac tions is the presenc e of an AGN in 
the galaxy NGC 7319 (jHuchra et al.lll982l) and the pres- 
ence of neutr al and molecular gas l ocated mainly outsid e 
the galaxies (jWilliams et al.l l2002l iLisenfeld etahl 120021) . 
UV observations with the GALEX satellite have provided 
a detailed picture of the unobscured component of recent 
star-forma tion in the group, including the tidal features 
(|Xu et al.l [2005). 

In the context of understanding the physical processes 
controlling star formation in groups in general and in SQ 
in particular the FIR spectral region contains crucial infor- 
mation. Specifically, IR emission from dust is fundamental 
for an understanding of the amplitude and distribution of 
star formation since it traces the luminosity of young mas- 
sive stars highly obscured by their parent molecular clouds 
at UV and optical wavelengths. The first infrared image 
of SQ, taken by the ISOCAM instrument on board the In- 
frared Space Observatory r evealed the pre sence of a star- 
burst, SQ A, in the IGM ()Xu et al.|[l999l ). possibly trig- 
gered by the o ngoing cohision between the intruder galaxy 
and the IG M (IXu et al.l l2003D. The later GALEX imaging 
of SQ (Xu et al.ll2005D also detected SQ A in the UV. Over- 
all, the distribution of UV emission in SQ measured by 
GALEX shows that most of the recent unobscured star for- 
mation has preferentially occured in the periphery regions 
or in the IGM. Further evidence for this is the discovery us- 
ing HST of widespread young star clusters distributed over 
the tid al debris and surrounding area by iGallagher et al.l 
(|2006( ). It is important to check whether this apparent 
shift from the main disk of the galaxies (wher e most star 
formation occurs in isolated galaxies, see e.g. iLeon et al.l 



l2008t) is also shown in the obscured component of star for- 
mation, and how the pattern of total star formation in the 
IGM is related to the morphological distribution of gas in 
the different temperature ranges. 

A further motivation for st udying SQ i n th e in- 
frared are theoretic al studies by iDwek etahl (|1990() and 
iMontier et al.l (|2004[ ) which have predicted that even small 
amounts of dust in the hot virialised IGM could pro- 
vide an important cooling mechanism via inelastic gas- 
grain collisions, with the radiation appearing in the FIR. 
This has prompted searches for a FIR counterpart to the 
X-ra y emitting iiitracluster mediurn in several rich clus- 
ters (IStickel etal.lll998l IStickel et al.]|2002l . iBai et al.ll2007l 



iKitavama et al. 20091 ) which however have thus far yielded 
no unambiguous detection. This may be attributed to the 
very low abundance of grains predicted on the basis of 
realist ic estimates of sour ces and sinks of grains in the 
IGM (jPopescu et al.|[20Q0l ). the expected similarity of the 
FIR colors of the collisionally heated dust emission compo- 
nent with photon-hea ted diffuse dust in foreground cirrus 
(jPopescu et al.l [20001 ) . and the problem of spati al confu- 
sion w ith star forming galaxies in the cluster (Q uillen et al.l 
Il999f ). SQ is a good object to search for this phenomena 
since its angular size is small enough for accurate pho- 
tometric measurements of extended emission and, at the 
same time, it is possible to separate emission from discrete 
star-forming sources, such as the constituent galaxies and 
objects like SQ A. Finally, the detection of dust FIR emis- 
sion, combined with gas mass measurements to obtain the 
dust to gas ratio, is a sensitive way of probing the metal- 
licity, and hence the origin, of the various gas components 
seen in the group in situations where optical nebular line 
diagnostic are weak or absent. 

The most extensive previous s tudy of the FIR emission 
from SQ was that by Xu et al... (2003) in the wavelengths 
range 11 to lOO/im using the Infrared Space Observatory. 
Apart from strong detections of the AGN galaxy NGC 
7319 and the foreground galaxy NGC 7320, significant de- 
tections were obtained in the FIR of the intergalactic star 
formation region SQ A and the source SQ B located on 
the tidal feature associated with NGC 7319. Intriguingly 
the ISO images also hinted at the presence of FIR emission 
associated with the X-ray emitting shock ridge, which was 
argued to be evidence for collisionally heated dust embed- 
ded in the hot shocked gas. 

In this paper we present deep imaging of the FIR emis- 
sion from SQ taken with the Spitzer Space Telescope. 
These data have superior angular resolution and sensitivity 
compared to the previous ISO study and moreover extend 
the wavelength coverage longwards to 160/J,m, thus pro- 
viding much deeper investigation of any cold dust compo- 
nents, embracing the expected spectral peak of the emis- 
sion from photon- and collisional-heated emissions. The 
data is used to characterize the pattern of star formation 
in SQ and to investigate possible infrared counterparts of 
the X-ray emitting structures. The Spitzer observations 
and data reduction are described in Sect. 2, where we also 
collate multiwavelength data needed for our analysis. In 
Sect. 3 we described the morphology of the dust emission 



Among the five galaxies that were observed for the first time by Edouard Stephan in 1877, NGC 7320 was later found to be a foreground 
dwarf galaxy. A nearby sixth galaxy, NGC 7320c, shows radial velocities compatible with SQ. Therefore the group is still a quintet even if the 
original definition has changed 
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of SQ in relation to the multiwavelength data set and also 
describe a novel fitting technique to the lower angular res- 
olution FIR maps which we use to separate the emission 
from the main emitting structures. After extracting pho- 
tometry from these structures across MIR/FIR range in 
Sect. 4, we fit the MIR/FIR SEDs with models in Sect. 5 
to elucidate the mechanisms powering the observed dust 
emission of the sources in SQ. This information is then 
used to quantify star formation rates in Sect. 6 where 
we also quantify corresponding gas reservoir available to 
fuel the star formation. In Sect. 7 we quantitatively dis- 
cuss the nature of star formation in and outside the con- 
stituent galaxies of SQ, together with the related question 
of the extent to which collisional heating plays a role in 
determing the thermodynamic properties of the IGM in 
SQ. In this paper, we assume a distance from the group 
equal to 94Mpc, corresponding to a systemic velocity of 
V — 6600 km/s and assuming Ho = 7Gkm/s/Mpc. 

2. OBSERVATIONS AND DATA REDUCTION 

In this work we have used Spitzer maps from two Guest 
Observer (GO) programs making use of continuum data at 
24, 7 and 160^m from the MIPS instrument ( Rieke et all 
|2004[) . From the GO #40142 (PI: Applcton, P.) we took 
the 24/um map and from the GO #3440 (PI: Xu, K.) the 
70 and 160/zm maps. In addition, a n image of SQ at 
8/Ltm, taken with the IRAC instrument (jFazio et al.ll2004l ) 
and downloaded from the Spitzer Science Center (SSG) 
archive, was also used. In the following we provide details 
of the data reduction for the Spitzer FIR maps. Th e 24/Ltm 
map has been already presented by iGluver et al.l (|2010( ) 
and we refer to that paper for technical details about the 
data preparation. 

2.1. MIPS 70/im and 160/im data reduction 

The MIPS 70/im and 160/zm data are from Spitzer 
pipeline version SI 1.0. 2. In the pipeline default mosaic 
images, there are a few residual instrumental artifacts 
that are particularly noticeable at 70/im. Because of this, 
these pipeline images were not used in this work. In their 
place, new images were used in which the instrumental ar- 
tifacts were reduced by performing additional data reduc- 
tion steps on the "basic calibrated data" (BCD) frames. 

Specifically, for the 70/im data, we observed an overall 
signal drift in time. To remove this, we masked out those 
pixels of each of the 468 non-stim BCD frames that are 
within a radius of 125" of the SQ center (RA=339.0181d, 
DEC=33.969183d; J2000). For each BCD frame, a median 
was calculated from all unmasked pixels. These medians 
were plotted as a function of the BCD frame index (1 to 
468). The resulting plot shows a clear discontinuity at 
frame index 313. We fit the sections prior to and post 
this discontinuity separately with a cubic spline function 
of order 1 or 2. The two fitted curves were connected to 
form one curve covering all the BCD frames. After be- 
ing normalized by its mean, this curve was divided into 
the index-ordered BCD frames to remove the signal drift 
in time. The next step was to create a sky flat image by 
median filtering only unmasked data for a given detector 

Throughout this paper all the S/im image s presented have had the stellar component of the emission subtracted using the relation 
(S^im, dust) = F^(8fj,ni) - 0.232F„ (3.6/im) lIHelou et al.|[200l) . The S^im Spitzer band also contains emission from rotational hydrogen 
lines which however we estimate in Sect. 5 to be unimportant in relation to the PAH and dust continuum emission. 



pixel. The resulting sky flat image was normalized by its 
mean, and subsequently divided into each unmasked BCD 
frames. Finally, we used Spitzer MOPEX tool to mosaic 
these improved BCD frames into our final image used in 
this paper. With a much flatter sky background, our own 
mosaic image is significantly better than the pipeline coun- 
terpart. 

For the MIPS 160/ini data, a similar procedure was used 
with all 522 non-stim BCD frames. In this case, the detec- 
tor signal drifts in time differ significantly among individ- 
ual readout modules. As a result, our signal drift removal 
was attempted on per readout module basis. 

Background subtraction has been performed by fitting 
a tilted plane to the maps, after masking a large area cov- 
ering the main group and, in the case of the 160/im map, 
the galaxy NGC 7320c (that lies outside the field of view 
shown in Fig. [T|). The total area used to estimate the 
background is about 50% of the entire maps at both wave- 
lengths. These areas are well outside the group emission, 
allowing good background measurements. 

2.2. Collected multiwavelength data 

In this paper we have also made use of a large set of 
multiwavelength data to derive physical quantities or for 
morphological comparison with the emission seen on the 
infrared maps. Specifically we utilized the SDSS maps 
of SQ, obtained from the SDSS Data Archive (data re- 
lease 7), the GALEX FU V map, which have been pre- 
sented in IXu et al.l (l2005l). the XMM- NEWTON soft X- 
ray map from Trinchieri et al.l (|2 005[ ). the VLA 21 cm 
line map from [Williams et al. (2002), the Ha maps from 
IXu et all (|1999(), the IRAC 3.6/im map, obtained from the 
SSC archive, and the CO radio observations presented in 
iLisenfeld et all (l200l . 

3. MORPHOLOGY OF DUST EMISSION IN SQ 

In Fig. |T] (lower panel) we show all the Spitzer broad 
bands (from IRAC and MIPS) where the signal is domi- 
nated by dust emissiorP^: 8/im, 24/im, 70/im and 160/im. 
These can be compared with the upper panel showing the 
SDSS r-band map, the GALEX FUV map, the VLA radio 
21 cm line map and the XMM NEWTON soft X-ray map. 
On each map, crosses identify the galaxy centres. 

The IR morphology of the galaxies seen in the IR maps 
is markedly different from the optical morphology. This is 
true not only for the early type galaxies but also for the 
late type galaxies which exhibit remarkably little infrared 
emission from the main bodies of the galaxies. The only 
exception is the emission from the foreground galaxy NGC 
7320 which is quite symmetrical in all the Spitzer bands, 
having a filled disk of emission commonly seen in local 
universe field galaxies, compatible with the optical/UV 
appearance not showing any sign of interactions (compat- 
ible with its not being a member of SQ). 

The infrared emission from NGC 7319 is dominated 
by an unresolved nuclear source, presumably from the 
Seyfert 2 AGN, that is particularly prominent at 24 and 
70/im. At 8/im one can also clearly see emission from the 
disk of the AGN host galaxy. As can be seen from the 
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MIR maps, dust emission from compact star formation 
regions are detected all over the group and especially on 
the elongated features of the intruder galaxy NGC 7318b 
which are most prominently delineated in the UV. Here 
the similarity is strongest between the UV and the 8/im 
band though the most prominent discrete sources are also 
clearly seen at 24 and 70/im. The most prominent such 
MIR/FIR source is SQ A, the star formation region lo- 
cated to the north of NGC 7318b, already detected by ISO, 
which can also be seen at 160/im. From optical spectra 
(IXu et al.' '2003 D and rad io observations (jLisenfeld et al.l 
12002^ Wilhams et al.ll2002D it is known that star formation 
in this region is associated with gas at radial velocities 
corresponding to both the intruder galaxy and the IGM of 
the group. Several further compact star formation regions 
are located on the southern arms of NGC 7318b. The 
brighter sources, HII SE and HII SW, are also detected 
on the 70/im map but not clearly seen on the lower reso- 
lution 160/im map. Two other bright MIR/FIR emitting 
regions detected on the Spitzer maps are SQ B, a star for- 
mati on region located on the optical "young " tidal tail 
(see ISulentic et all l200lD . and a source, HII N, located 
about 40" towards the north of the AGN galaxy. 

There is no clear morphological counterpart in the in- 
frared to the shock region, defined here by the ridge of 
emission that can be seen on the X-ray map. Neverthe- 
less the 70/im map and more particularly the 160/im map 
show enhanced emission towards the peak of the X-ray 
emission. The ratio of the MIR 8 and 24/im emission to 
the 160/im emission at the same position appears rather 
low compared to the discrete sources associated with star 
formation such as SQ A. A previously undetected feature is 
an extended FIR emission component, spatially coincident 
with the main part o f the group X-ray halo (as defined in 
iTrinchieri eralll2005D . The appearance of the correspond- 
ing MIR emission on the higher resolution 8 and 24/im im- 
ages suggests that the extended FIR emission may at least 
in part be clumpy rather than uniform, possibly indicating 
the presence of faint star formation regions far away from 
the centers of the galaxies. 

In order to quantify the morphology and brightness of 
the FIR emission seen towards the X-ray emitting halo 
and shock regions, it is necessary to subtract from the FIR 
maps the most prominent discrete sources associated with 
star formation regions and galaxies. To do this we devised 
a FIR source fitting technique which we describe in sub- 
section [311] This source fitting technique also serves to fix 
the photometry and the extent of discrete sources, infor- 
mation that cannot be directly extracted from the maps 
due to the unknown level of mutual confusion. 

3.1. The FIR map fitting technique 

The fitting technique models a preselected set of the 
most prominent discrete sources as a sum of elliptical gaus- 
sian convolved with the PSF. Seven parameters are calcu- 
lated for each source: amplitude, the peak coordinates, 
the two gaussian widths, the axis rotation angle and the 
local background (included to avoid removal of any diffuse 
emission components). The fit is performed simultane- 
ously for the 10 brightest sources seen on the 70/tm maps, 
(see Fig. [5]): five compact sources (SQ A, HII SE, HII 
SW, SQ B, HII N), two sources to fit the emission from 



the AGN galaxy NGC 7319 (one for the central emission 
and one for a peripherical star formation region visible af- 
ter the removal of the first component), two for the fit of 
the foreground galaxy NGC 7320 (one for the fit of the 
diffuse emission and one for a compact source) and one to 
fit the emission peaked in the middle of the shock region. 
The fit to the 70/im map is performed first, keeping all the 
fitting parameters as free variables. This is followed by a 
constrained fit to the 160/im map in which the relative 
position of the sources are fixed to the values obtained 
at 70/im. In addition, a further constraint was that we 
kept the same shape and axis orientation for the five com- 
pact sources, as inferred by the 70/im fit, allowing only 
a size change (to take into account the potentially more 
extended distribution of cold dust emission). This strat- 
egy was adopted because the higher resolution 70/im map 
places the strongest constraints to the position and mor- 
phology of sources in the FIR. We did not use the highest 
resolution 8/im map for this purpose since, whereas the 
emission detected at 70 and 160/im is produced by the 
same kind of solid dust grains, the emission at 8 /im is 
dominated by line emission from PAH molecules. Model 
images from the fitting technique are shown in Figs. [5] and 
[21 The best fit parameters and the inferred total fiux den- 
sities are shown in Table [TJ At both 70/im and 160/im the 
best fit model images are remarkably similar to the origi- 
nal maps. On the "deconvolved" maps in each figure, it is 
possible to see the contribution from each gaussian to the 
final map. Interestingly enough, the source at the position 
of the shock is much more predominant at 160/im than at 
70/im, confirming the original impression that the emis- 
sion in the shock region is brighter at 160/im. It is also 
noteworthy that the position angle of the model source 
at 160/im is aligned with the north-south orientation of 
the X-ray emitting ridge whereas at 70/im no such align- 
ment is apparent. At both 70 and 160/im the east-west 
width of the fitted elliptical gaussians at the shock po- 
sition {FWHM w 60") is larger than that of the X-ray 
shock ridge {FWHM k, 20"). This indicates that the in- 
tegrated emission is not necessarily entirely composed of 
emission from the shock ridge. A full description of the 
source fitting technique is given in Appendix 1X1 

3.2. The FIR residual maps 

To understand how the emission peaked in the shock 
region and the extended emission are distributed on the 
maps, we created FIR residual maps where the emission 
from all the sources fitted by the FIR map fitting tech- 
nique, with the exception of the source associated with the 
shock, have been subtracted. The contours of these FIR 
residual maps are shown in Fig. SI overlaid on HI, X-ray 
and FUV maps. As one can see, the emission on the FIR 
residual maps is uncorrelated and perhaps even anticor- 
related with the HI distribution but well correlated with 
the soft X-ray flux. The FIR emission, as already seen 
on the original maps, peaks in the middle of the shock 
region and its overall extent is similar to the X-ray halo 
emission (similar to what wa s seen in the low e xcitation 
pure rotational lines of H2 bv lCluver et al.|[2"010[ ). At first 
sight this finding supports the idea that collisional heating 
is producing the observed FIR emission. However a large 
part of the residual FIR emission covers areas emitting 
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significant luminosity at UV wavelengths. The presence 
of these radiation sources complicates the interpretation 
of the dust emission in the shock region as well as for the 
extended emission (see Sectl?!). 

4. SQ SPITZER MAP PHOTOMETRY 

Precise photometry of all these different emitting regions 
in SQ is required in order to elucidate the physical mech- 
anisms that power dust emission and the related scientific 
implications. In the FIR the photometry is derived from 
the source fitting procedure (for the discrete sources) and 
from the residual maps (for the shock and the extended 
emission components). In the MIR the high resolution of 
the MIR 8 and 24/im maps allows a straightforward ex- 
traction of fluxes of the corresponding regions in SQ using 
aperture photometry. All results are summarized in Table 

m 

4.1. Star formation regions and galaxies 

The FIR map fitting technique, described previously, 
allows a precise measure of the fiux coming from the com- 
pact star formation regions, the AGN galaxy NGC 7319 
and the foreground galaxy NGC 7320, all well modelled by 
convolved elliptical gaussians (although some sources are 
not detected at 160/im). This technique allowed us not 
only to obtain the total source fluxes but also to derive 
accurately the source extent. In Fig. [51 as an example, we 
show the contours of the "deconvolved" 70/im and 160/im 
emission at the position of SQ A overlaid on the 8/xm and 
24/im maps. The FIR emission has a MIR counterpart 
that peaks in the areas where FIR is higher. This is gen- 
erally true for all the fitted FIR sources and it validates 
the use of apertures for the photometry of the compact 
sources at MIR 8 and 24/im whose sizes are derived by 
the FIR fitting technique. Specifically, at 8 and 24/im we 
used elliptical apertures having the same axial ratio and 
orientation as the 70/xm best fit elliptical gaussian axis 
and semi- axis lengths equal to 2.17(T7o^m (this area in- 
cludes 90% of an elliptical gaussian total flux). We took 
sizes based on the 70/<m map fit because all the compact 
sources are clearly seen on that map. These apertures are 
depicted in red in Figs. in]and[71 The local background for 
the aperture photometry has been estimated on regions 
located nearby the central source, where other peaks of 
emission are not clearly seen (see green circles on Figs. [5] 
and [7]) . In this way we are confident that the background 
level has not been overestimated due to contamination by 
surrounding sources. We didn't apply aperture correction 
for the photometry at 8/<m because this is close to unity for 
the chosen apertures (typical aperture size about 10 pix- 
els), as reported by the IRAC data Handbook. At 24/im 
the chosen apertures typically delimit the first outer ring 
of the PSF. In this case we applied an aperture correction 
equal to 1.16 (MIPS data Handbook). For the MIR pho- 
tometry of the galaxies NGC 7319 and NGC 7320 we used 
large apertures covering most of the emission from these 
objects and we didn't apply aperture corrections because 
of the large integration area. 

The uncertainties on the MIR photometry are given by 
the quadratic sum of the following contributions: 1) error 

The presence of a compact X-ray source (see Fig. [TJ and a con 
luminosity AGN. 



on the aperture correction; 2) flux calibration uncertainty; 
3) background fluctuations. Since we used elliptical aper- 
tures, instead of typical circular apertures, we rather con- 
servatively assumed that the relative error on the aperture 
correction is 10% (note that this error is applied only to 
aperture photometry at 24/im). The flux calibration rel- 
ative uncertainty is equal to 4% (IRAC and MIPS Hand- 
book) while the error introduced by background fluctua- 
tions is derived from the variance of the background mean 
values in the several areas we selected nearby the sources. 
The measured fluxes and uncertainties for both MIR and 
FIR wavelengths are given in Table [51 

Since we are using different methods to derive the source 
fluxes on the MIR maps and the FIR maps, it is impor- 
tant to check that these two photometric techniques give 
consistent results. To verify this, we performed the source 
fitting technique on convolved MIR maps and derived the 
source fluxes exactly as we did for the FIR maps. Details 
of this test are shown in Appendix \^ In Fig. [8] we plot 
the ratio of the 24/im fluxes inferred by the fltting tech- 
nique and by aperture photometry for each source. For all 
the sources, except for HII SW, the difference between the 
fluxes obtained using different methods is less than 20%. 
The discrepancy between the fluxes for HII SW arises be- 
cause of the presence of a compact 24/im source at the 
nearby nucleus of NGC 731843- This source is included in 
footprint of the gaussian returned by the fitting procedure 
but is outside the integration area used for the aperture 
photometry. 

4.2. Shock region 

Fig. [31 show the east-west profiles, along a line passing 
through the center of the shock ridge, extracted from the 
FIR residual maps, the X-ray and the VLA 21 cm radio 
continuum maps convolved to the resolution of the 160/im 
map {FWHM w 40"). From these profiles, one can see 
that the FIR width is significantly larger than the shock 
ridge width as seen both in the X-ray and Radio. Using 
the current data it is impossible to say if this discrepancy is 
due to confusion with fainter unrelated infrared sources or 
to a systematic change in the width of the emitting region 
of the shock between the X-ray and the FIR. Nevertheless 
we estimated the flux coming from the shock ridge by fit- 
ting the FIR residual maps with a simple two component 
model: a PSF convolved uniform ridge, used to fit the dust 
emission in the shock region and whose size (20" x 80") 
were derived from X-ray data, plus a uniform component. 
In this procedure the level of the fitted uniform compo- 
nent is influenced by the fact that the actual FIR source is 
more extended than the X-ray source. Taking into account 
the ambiguity in identifying all the flux from the more 
extended FIR emission with the shock, we have assigned 
extremely conservative flux uncertainties. Specifically, the 
upper and lower limits defined by the quoted error bars in 
Table [21 are the fluxes contributed by the uniform emission 
component underlying the solid angle of the shock ridge 
(convolved with the PSF) in the cases the uniform emis- 
sion component has a brightness equal to twice or zero 
times the value given by the fit. 

To measure the MIR 8 and 24/im emission from the 

radio source (see IXu et al.ll2003l) suggest that this source is a low 
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shock region, we simply integrated the emission in the 
same rectangular area used before to define the shock 
ridge component in the fit of the FIR residual maps. The 
adopted apertures are shown in Figs. [SlandlTl Before the 
integration we masked all the areas inside the apertures 
we used for the photometry of other sources, because the 
MIR emission in that areas is mainly connected with the 
corresponding sources that we fitted and subtracted from 
the FIR maps before the measure of the shock region flux. 
For the fraction of the masked regions that falls inside the 
shock region integration area we assumed that the surface 
brightness is equal to the average brightness on the other 
regions inside the aperture. The background level has been 
measured on areas around the rectangular aperture where 
no peaks of emission are clearly seen. The quoted error on 
the integrated fluxes is the sum of the contributions due 
to background fluctuation and flux calibration error. 

4.3. Extended FIR emission 

We have measured the amount of extended flux on the 
FIR residual maps within a radius of 90" from the shock 
center, an area approximately equ al to that covered by th e 
X-ray HALO region as defined in iTrinchieri et all (|2005[ ). 
We performed this flux measurement in the following way. 
First we constructed radial curves of growth of the inte- 
grated emission on the FIR residual maps starting from 
the shock center. These curves of growth are shown in 
Fig. [TUl As one can see, the integrated emission continues 
to grow somewhat beyond the 90" radius. However we 
considered only the flux within this limit because it can 
be directly related to the X-ray Halo extent. We took the 
curve of growth values for the integrated fluxes at 90" and 
we subtracted the contribution from the shock region, es- 
timated in Sect. 14.21 Similarly as before, the uncertainties 
on the fluxes are mainly due to the mutual contamination 
between shock ridge and extended emission. Therefore we 
assumed the same conservative errors that we assigned to 
the shock region FIR fluxes. For the estimate of the ex- 
tended MIR emission we have used an analogous method. 
We constructed radial curves of growth, starting from the 
shock region center, after having masked all the compact 
sources and galaxies whose photometry has been described 
in Sect. l4.1l fnote that the shock region is not masked). For 
the calculation of the MIR curves of growth we took into 
account the missing areas, those that we masked, assum- 
ing that their brightness is equal to the average brightness 
inside the circular annuli passing through them. The de- 
rived curves are shown in Fig. (TU] Exactly as for the FIR 
measurements, we took the value of the curve of growth 
at 90" and subtracted the emission from the shock region 
in order to obtain the integrated emission from the ex- 
tended area corresponding to the X-ray halo. The flux 
uncertainties, in this case, are derived summing quadrat- 
ically background fluctuation, calibration errors and the 
error on the shock ridge flux. 

5. MODELLING THE INFRARED SEDS 

The spectral energy distribution (SED) of dust emission 
is determined by the heating mechanism, the intensity and 
color of the radiation fields (for the case of photon heat- 
ing), the temperature and density of the hot plasma (for 
the case of coUisional heating) and by the amount, size dis- 



tribution and chemical composition of the emitting grains. 
Ideally one would solve for the distribution of photon 
sources, hot plasma and grains using a self consistent radi- 
ation transfer analysis to fit the entire X-ray/UV/Optical- 
MIR/FIR SED for each source, analogous to the treat- 
ment of pho ton heated dust in disk galaxies and starburst 
galaxies fsee lPopescu fc Tuffs I ([2010) for a recent review). 
However, due to the extra dimension of collisional heating 
and the unusual geometry of the optical/UV and infrared 
emission from SQ, this approach would require a very in- 
dividualized treatment which is well beyond the scope of 
this paper. We therefore adopt here a hybrid approach 
fitting just the MIR/FIR SED with superpositions of dust 
emission templates where each template is appropriate for 
specific dust emission regions: Hll/photodissociation re- 
gions (FDR), diffuse photon-powered dust emission, AGN 
torus emission, collisionally heated dust embedded in hot 
X-ray plasma. Each of these templates is calculated self 
consistently in terms of physical input parameters, as de- 
scribed in detail in Appendix |B] In this approach the am- 
pl itude of the HII/PD R te mplate, derive d from the model 
of lDopita etHI (|2005D and lGroves et al.l (^2008) and shown 
in Fig. [TT](left panel), quantifies the obscured component 
of on-going star formation, whereas the amplitude of the 
diffuse dust emission template (shown in Fig. I12p quan- 
tifies dust emission powered by longer range photons to 
which older stellar populations can also contribute. In 
the latter case a family of SEDs is calculated according 
to the strength and color of the radiation field which are 
both expected to vary according to radiation transfer ef- 
fects and the relative contribution of young and old stel- 
lar populations. A family of SEDs is also calculated for 
the collisionally-heated dust emission template (Fig. [T51) 
corresponding to a range of plasma parameters and to 
different grain size distributions thus accounting for the 
expected effect of the plasma on the grain size distribu- 
tion. For both the photon heated and collisionally heated 
dust emission templates the stochastic emission from im- 
pulsively heated grains, which is important to determine 
the MIR emission, is calculated. In the case of the AGN 
template SED, we have u sed the existing self consistent 
model of iFritz et all ()2006[ ) (right panel of Fig. HI]). 

The fits were done by minimizing for each source, 
defined as: 

/ TT'modcl i^obs \ ^ 

where F°^^ and ai respectively are the observed fiux den- 
sities and their associated uncertainties , and F™°'^°^ is the 
corresponding model prediction which is related to a given 
theoretical spectra through i^'modoi _ KF^. The color 
correction K for each Spitzer band is calculated accord- 
ing to the formulas provided by the IRAC and MIPS data 
handbooks (see Appendix |D]). It is particularly impor- 
tant for the 8/im band where the emission spectra varies 
quickly along the bandwidth because of PAH line emis- 
sion. The color correction can be up to 80% at 8fj.m, while 
it is generally less than 10% for the other bands. 

5.1. Star formation regions 

In this section, we describe the SED fits for the star 
formation regions SQ A, HII SE, HII SW and SQ B. We 
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fitted all the observed SEDs as the superposition of two 
components. The first component is the PDR/HII region 
dust emission template, that helps to fit the emission from 
warm dust located very close to young stars. The sec- 
ond component is the diffuse photon heated dust emission 
template that fits the emission from diffuse dust near star 
formation regions. This diffuse dust is heated by a combi- 
nation of UV photons escaping from the PDR/HII region 
and any ambient large scale UV/optical radiation field per- 
vading the region. We performed the fit varying four free 
parameters: the amplitude of the HII region/PDR SED 
template (determined by the parameter F^^, the fraction 
of 24/im flux contributed by this component); the diffuse 
dust mass Mdust ; the strength and the color of the diffuse 
radiation field, determined by two parameters, Xisrf and 
Xcoior (see Appendix [B] for details). As one can notice, 
there are not degrees of freedom in these fits. Therefore 
it is not possible to estimate the goodness of the fit (in 
the sense of the fidelity of the model) from a chi-square 
test. Nonetheless we estimated the error on the best fit 
dust masses and total dust luminosities from a multidi- 
mensional analysis of near to its minimum. The uncer- 
tainty on the dust mass is the minimum dust mass vari- 
ation that gives Ax^ = — Xmin values always greater 
than one, independently from the values of all the other 
parameters. The total luminosity error bar is determined 
by the lowest and highest values of the total luminosity in 
the subspace of fitting parameters determined by Ax^ < 1. 

The SED fits are shown in Fig. [H (note that for HII 
SE and SQ B we used the I60/im fluxes within the 70/im 
source size to limit contamination from neighbours). Ta- 
ble [3] shows the best fit parameters, together with the 
estimate of the uncertainty on the dust mass, the total 
infrared dust luminosity and the contribution to this to- 
tal infrared luminosity from optical heating of diffuse dust 
by the diffuse radiation field, UV heating of diffuse dust 
by the diffuse radiation field and localized UV heating of 
dust in PDR/HII regions. The predicted radiation fields 
needed to account for diffuse emission component are at 
least as strong as the local ISRF in the Milky Way with 
a strong variation in color. From the plots one can see 
that the 24/im emission is dominated by HII region/PDR 
emission, as expected, while 8/im and FIR points are gen- 
erally dominated by diffuse emission, in accordance with 
fully self-consistent r adiative transfer models of galaxies 
(jPopescu et al.|[2010l submitted). This is consistent with 
a picture where MIR emission is produced by regions very 
close to young stars whilst FIR emission and PAHs emis- 
sion come from dust illuminated by dilute rad iation fields, 
as observed in galaxy disks (e.g. iBendo et al.|[2008i) 

5.2. AGN galaxy NGC 1319 

The observed SED of the Seyfert 2 galaxy NGC 7319 
differs markedly from the SEDs we obtained for the other 
sources in SQ, for which the power emitted at 24/im is 
typically a factor ten or more smaller than that observed 
at 70 and 160/im. For NGC 7319 the amount of energy 
emitted at MIR wavelengths is comparable to the FIR lu- 
minosity. We first tried to fit the SED of this galaxy using 
the same HII region/PDR plus diffuse emission compo- 

Comparing the surface brightnesses in the — 05'{4) and 
("Main Shock" region), with the 8^m surface brightness derived from the data 
of the Spitzer S/xm band by II2 hne emission to be 10%. 



nents as we have done previously for the star formation 
regions (we didn't include a component of synchrotron 
emission, originating from the AGN, because an extrapola- 
ti on of its infrared lu minosity from the radio measurements 
of lAoki et al] (IT991 gives values which are six orders of 
magnitudes lower than the observed infrared luminosities). 
The best fit is shown in the left panel of Fig. [TS] The fit 
tries to reproduce the high 24/im fiux with the HII re- 
gion template but, doing so, completely overestimates the 
70/xm flux. A second attempt has been performed substi- 
tuting the HII regioii/PDR template with an AGN torus 
template (iFritz et al.ll2006t see Appendix iBt . As one can 
see from the right panel of Fig. [TSl the observed data are 
reproduced much better in this case. Of course the lim- 
ited number of data points does not allow any inference 
of the physical properties of the dusty torus or provide 
precise information on the radiation held. However the fit 
does provide an indication that the majority of the 24/im 
emission in this galaxy is powered directly by the AGN 
power law radiation. The fitted parameters and the total 
dust luminosity derived from the AGN torus plus diffuse 
emission fit are shown in Table [3] 

5.3. Shock region emission 

The dust emission SED of the shock region contains MIR 
and FIR emission components which indicate an inhomo- 
geneous structure of the emitting region. As one can see 
from the middle panels of Fig. HI the emission on the 
FIR residual maps shows a rough correlation with X-ray 
emission, especially near the center of the shock ridge, sug- 
gesting, as mentioned before in Sect. 13. 2[ that collisional 
heating of dust embedded in the hot (T w 3 x 10^ K) 
shocked plasma powers at least some of the observed emis- 
sion. However the presence of significant 8/im emission 
rules out the possibility that the observed dust emission 
is purely coUisionally heated because PAH molecules are 
not expected to survive in a medi um shocked by a wav e 
travelling at more than 125 km/s ()Micelotta et al.ll2010[ ). 
The PAH emission must therefore either be physically un- 
related to the shock, coming from another region along the 
line of sight, or, if associated with the shock, belonging to 
colder gas phases embedded in the X-ray plasma. The lat- 
ter scenario woul d be consistent with the detection of H2 
molecular lin£f| (Apple ton et"alll20Q6l . ICluver et al.ll2010D 
and Ha (|Xu et al. 1999[) emission from the shock region. 
Peculiar ratios between the PAH line features have been 
detected towards the shock, which suggests an enhanced 
fraction of neutral a nd large PAHs compa red to typical 
galaxy environments (' Guillard et al.l (j2010 ^). 

The multiphase nature of the gas in the shock region im- 
plies that there are at least four potential sources of dust 
emission in this region: 1) diffuse dust coUisionally heated 
in hot plasma; 2) dust in a colder medium and heated by a 
diffuse radiation field; 3) HII/PDR dust emission from op- 
tically thick clouds with embedded star formation regions; 
4) cold dust emission from optically thick clouds without 
embedded star formation. 

In practice, it is not possible to distinguish model pre- 
diction for 2) and 4) over the wavelength range of the 
currently available data, since the only difference would 

05'(5) lines, derived from the fluxes given in Table 2 of lCluver et~all lIMTol ') 
iven in Table [2] of this paper, we estimate the contamination 
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correspond to a very cold emission component from the 
interior, self-shielded regions of the optically thick clouds, 
which will only become apparent at longer submm wave- 
lengths. Therefore, we only consider explicitly here com- 
ponents 1) - 3), ignoring component 4). Furthermore, the 
four available data points are insufficient to simultaneously 
fit these three potential sources of dust emission, especially 
the diffuse photon heated and collisionally heated com- 
ponents, which are both predicted to peak at FIR wave- 
lengths. Therefore, we decided to perform the SED fit in 
two different ways corresponding to the two opposite cases, 
where collisional heating is either responsible for the entire 
FIR emission or is completely negligible. We followed this 
approach in order to understand which mechanism is pre- 
dominant in powering the diffuse dust emission. First, wc 
fitted the shock region SED as a superposition of two com- 
ponents: the HII region/FDR template, to fit MIR emis- 
sion possibly associated with star formation regions and 
a collisionally heated dust SED template (see Appendix 
IB. 0.51) . This model is appropriate for the case where pho- 
ton heated diffuse FIR emission is negligible compared to 
the collisionally heated dust emission. Given the plasma 
physical properties (fixed by the X-ray emission, as de- 
scribed in Appendix [C]), the free parameters in this fit are: 
the amplitude of the HII region template and the colli- 
sionally heated dust mass. We performed these fits for 
the highest and the lowest densities admitted by the X- 
ray data, however, the final results are quite similar. In 
the upper panels of Fig. (THlwe show the best fit obtained 
for the highest considered particle density n ~ 0.016 cm^"^ 
and temperature T = 3 x 10^ K. As one can see, the col- 
lisionally heated component (dashed curve) is sufficient to 
reproduce the observed data except the 8/xm point which, 
as already remarked, is too high to be due to PDR/HII re- 
gions, requiring a diffuse photon-heated contribution. The 
model curve shown (fitted in this case just to the three 
longer wavelength data points) was calculated assuming a 
power law grain size distribution with exponent k — 2.5, 
the expected value in the case of equilibrium between dust 
injection (with a st andard fc = 3.5 interstellar distribution; 
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19771) and sputtering of grains in hot plasmas 
19901) . The dust mass required to produce the 



observed FIR fiux densities for the case where collisional 
heating, rather than diffuse photon heating, is considered 
is 2 X 10^ Mq. 

The second fit was performed using the HII region/FDR 
template and the photon- heated diffuse dust component. 
This combination of SED templates represents the case 
where the entire emission is powered by photons and col- 
lisional heating is negligible. Using this combination of 
SEDs, one can fit all the data points, including the S/um 
fiux. The best fit curve is shown in the lower panels of 
Fig. 1161 In this fit, the diffuse dust component dominates 
the emission in all the Spitzcr bands, even at 24//m, where 
the contribution of PDR/HII regions is normally predom- 
inant. The fitted parameters are given in Table |31 We 
note that the diffuse radiation fields needed for the fit are 
colder than those needed to fit the star formation regions 
in Sect. O 

In conclusion, the FIR part of the SED from the shock 
can be fitted by a cold continuum component either from 
collisional heating, which is cold due to the limited ambi- 



ent density of plasma particles and the underabundance of 
small stochastically heated grains, or from photon heating, 
which is cold due to the low ambient density of photons 
coupled with a low UV to optical ratio. The only real 
discriminant favouring photon heating would be the S^m 
emission which, as a PAH tracer, cannot be explained by 
collisional heating. However it is not entirely clear, on the 
basis of the current data, that the 8/im emission really 
does all originate from the shock region. 

5.4. Extended emission 

In a completely analogous way as for the shock region 
emission, we tried to model the observed extended emis- 
sion SED by including a collisionally heated dust emis- 
sion component. This is motivated by the fact that the 
FIR residual maps show the presence of extended emis- 
sion in the area of the brightest compon ent of the X-ray 
halo emission (the so called "HALO" in iTrinchieri et al.l 
120051 ). In the upper panels of Fig. [T7] we show the SED 
fit performed adding two components: HII regions plus 
collisionally heated dust emission. In this fit the adopted 
plasma physical parameters, derived in Appendix [C] are 
n — 0.001 cm~'^ and T = 6 x 10^ K and we again as- 
sumed k — 2.5. Using this combination of SEDs, one can 
roughly reproduce the SED longward of 24/zm, but the 
collisionally heated component dominates only the 160/Ltm 
fiux. The dust mass inferred from this fit (constrained, as 
in the case for the model fit to the shock SED incorporat- 
ing collisional heating, only by the three longer wavelength 
points) is 1.0 X 10® Mq. 

In the lower panel of Fig. [T71 we show the fit performed 
adding HII region and diffuse photon heated emission. As 
before for the shock region, using these two components 
one can fit the entire spectra including the 8/im point. In 
this fit, HII region emission dominates at 24 and 70/xm 
while the diffuse emission is responsible for the 8 and 
160/im fiuxes. The fitted parameters are shown in Table 
[3l As for the shock emission, we note that the extended 
emission component requires a very cold FIR component 
from the fit even though, observationally speaking, this is 
a somewhat less robust conclusion, since it really needs to 
be confirmed with longer wavelength photometry. 

5.5. Foreground galaxy NGC 7320 

For the SED fit of the foreground galaxy NGC 7320, we 
used again the usual combination of HII region plus dif- 
fuse photon heated dust emission components. The best 
fit is shown in the left panel of Fig. [TSl The relative 
contribution of HII regions to the 24/im fiux is reduced 
compared to the star formation regions, whose SED fits 
have been shown before. This might be physically under- 
standable since, in this case, we are fitting the emission of 
a galaxy as a whole, including all the diffuse interstellar 
medium containing small dust particles whose stochasti- 
cally heated emission ca n account for a majo r part of the 
total MIR emission (e.g. iPopescu et aLll2000t ). 

5.6. Infrared source north of NGC 7319 (HII-N) 

We attempted to fit the observed emission for this source 
with a combination of HII region plus diffuse photon- 
heated dust emission. However, this source shows a very 
peculiar MIR to FIR ratio and, as a consequence, we have 
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not managed to fit the entire spectra. As shown in the 
right panel of Fig. [181 the FIR points are highly under- 
estimated by the best fit curve. If no systematic errors 
are present in our measurements, it might be that HII 
N is a distant background source and the intrinsic colors 
are heavily redshifted. This possibility is strongly sup- 
ported by the presence of a red spiral galaxy at the po- 
sition of HII-N, clearly visible on recently released high 
resolution HST color maps of SQ (Hubble Servicing Mis- 
sion 4 Early Release Observations, observers: K.Noll et al., 
available at http : / /hubblesite . org / newscenter / archive /( re- 
leases/2009/25/image/x/). 

6. QUANTIFYING STAR FORMATION RATES AND GAS 
MASSES IN SQ 

6.1. Star formation rates 

The measurements of UV, recombination line and dust 
emission from sources in SQ can in principle be used to de- 
rive star formation rates (SFR), provided proper account 
is taken of the absorption of the UV/optical photons by 
dust and subsequent re-emission in the MIR/FIR spec- 
tral regimes. Several authors have provided empirically- 
based re lations achievi n g this for spiral galaxies on scales 
of kpc. ICalzetti et al.l (j2007[ ) presented an Ha luminos- 
ity based star formation rate relation: SFR{MQyr~^) = 
5.3 X 10-42 [L{Ha)ohs + (0.031 ± 0.006)L(24/xm)]. In this 
relation the 24/im luminosity is used to esti mate the Ha 
luminosity obscured by dust. .Bigiel et al.l pOOS) com- 
bined GALEX FUV and Spitzer 24/Lim fluxes to obtain 
SFR per unit area maps of a sample of spiral galaxies. 
Like ICalzetti et al.l (|200lD . they used a UV SFR calibra- 
tion and they used the 24/im flux to measure the obscured 
UV flux. 

Simple application of the Calzetti relation to SQ is ham- 
pered by the lack of pure measurement of Ha emission 
driven by sta r formation. As one can see from Figs. 4 and 
5 from lSulentic et al.l ()200ll ). the shock region Ha emis- 
sion shows a diffuse component that reflects the north- 
south ridge seen in the soft X-ray regime. This diffuse Ha 
emission cannot be conside red for SF R measurements be- 
cause, as demonstrated bv IXu et al.l ([2003) from spectral 
line analysis, the Ha emission in the shock region is dom- 
inated by shock-excitation rather than star formation. In 
addition, it is not clear that the semi-empirical relations 
derived for spiral galaxies should apply to the sources in 
SQ. 

Therefore, we adopted a new approach, using only obser- 
vational indicators of star formation activity available for 
all sources in SQ, and utilizing the results of the fits to the 
dust emission SEDs given in Sect. jS) Our method to esti- 
mate SFRs is based on a UV - SFR calibratioqYl. Spec if- 
ically we adopted the calibration bv lSalim et al.l (|2007t ): 

SFRuY = 1.08 X IQ-^^Fuv Mg/yr (2) 
where i^uv is the FUV luminosity density in units of 
erg/s/Hz which would be observed in the GALEX FUV 
band in the absence of dust. Fjjy can be written as: 

T71 Indirect I TTiabs, local . rj^abs, diffuse 

where -F'uv°'^* is the directly observed unabsorbed com- 
ponent of FUV luminosity density, p^^'^°'^'^^ jg tj^g PUV 

UV emission traces star formation on time scales of ~ lOOMyr and, 
formation on smaller time scales. 



luminosity density locally absorbed by dust in star forma- 
tion regions and ^^^^''1'*^^'''= jg the FUV luminosity density 
absorbed by dust in the diffuse medium surrounding the 
star formation regions. 

We measured -Fu™'^' fo'" each source by performing aper- 
ture photometry on the GALEX FUV map in a completely 
analogous way to the photometry we performed on the 
Spitzer MIR maps, including the construction of the curve 
of growth after masking the galaxies NGC 7319 and NGC 
7320 and the star formation regions (see Fig. jTU]). The 
final flux densities, shown in col. 4 of Tabled include the 
correc tion for Galactic fore ground extinction {E{B — V) = 
0.079. iSchlegel et al.l [19981 A(FUV) = 8.24 * E(B - V), 
iWvder et al.ll2007D . Values for the obscured emission com- 



, 7-^abs, local j lyiabs, diffuse , j. j r j.i 

ponents i'uy and i'uy were extracted from the 

fits to the dust emission SEDs by noting that the total 
infrared luminosity emitted by dust and powered by UV 
photons can be written as: 
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ing, and 2^^t)S''i'ff"''<= jg the part of the diffuse dust lumi- 
nosity powered by UV photons, respectively tabulated in 
cols. 6 and 8 of Table [3] Since the intrinsic SEDs of the 
young stellar populatio n are rather flat at UV wavelengths 
(jKennicutt et al.lll998l) . it then follows that: 

r abs, local 



is the luminosity of dust emission in 
regions, dominated by UV photon heat- 



and 



F 



rTiabs, local 
UV 
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UV 
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abs, diffuse 
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where Ai' « 1.8 x 10^^ Hz is the UV frequency width. The 

nabs. local , 

'uv' + 



total obscured UV luminosity density F^y — F^ 

T-iabs, diffuse 



UV 



and SFRs are shown in cols. 6 and 7 of Table 

As a check of the consistency between our method to 
derive st ar formation rates an d the Ha 24/im SFR re- 
lation of lCalzetti et al.l (|2007f) , we used the latter to derive 
SFR for the compact star formation regions in SQ which 
are the objects closest resembling galaxies. The Ha fluxes 
have bee n measur e d on the interference filter maps pub- 
lished in lXu et al.l (|l999 t) while the 24/im fluxes are those 
derived by aperture photometry (Sect. 14. ip . The results 
are shown in cols. 2 and 3 of Table [H As one can see com- 
paring cols. 3 and 7 of that table, the SFR inferred with 
our method are consistent with the results found using the 
Calzetti relation (except for HII SW but in that case the 
SED fitting was performed without varying Xcoior because 
of the non detection of the source at 160/im.) 

6.2. Gas Masses 

In the absence of a homogeneous set of observations of 
gas tracers for the gas phases of interest, namely HI, H2 
and X-ray emitting plasma, we proceeded differently for 
the several sources in determining gas masses. For the 
star formation regions SQ A, HII SE, HII SW and SQ B 
we measured gas mass column densities from the HI map 

therefore, UV derived SFRs could not be appropriate for burst of star 
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published by IWilliams et alj (|2002D and the CO maps of 
iLisenfeld et all (|2002[ ). Specifically we measured the ayer- 
age atomic and molecular hydrogen column density in the 
obseryed areas close to the position of the starburst re- 
gions. The inferred gas mass column densities are shown 
in col. 9 and 10 of Table 0] (note that for some yalues 
only upper limits are ayailable and no CO obseryation are 
ayailable for HII SE). 

Gas in the shock region is mainly in the form of X-ra y 
emitting plasma and molecular gas (iGuillard et al.l 20091) . 
From the X-ray luminosities, measured by lTrinchieri et al.l 
(|2003), we determine the X-ray gas mass Mx « lO^M© 
(see Appendix [C] for details). Diyiding the total gas mass 
by the physical area coyered by the shock (330kpc^), we 
obtained the hot gas mass surface density shown in col. 
11 of Table 01 To determine the molecular gas mass 
surface density, we measured the ayerage of several ob- 
seryed position s withi n the shock area on the CO maps by 
ILisenfeld et all (|2002D to obtain ^(iJa) =9 Mq/pc^. As 
one can see from Fig. 1 of that paper, most of the ob- 
seryatio ns were perform ed in the upper part of the shock 
region. IGuillard et al] pOlO.) and ,Applcton ct al. (2006) 
found a cold and warm molecular gas surface mass den- 
sities in the central parts of the shock region respectively 
equal to 5 Mo/pc^ and 3.2 Mo/pc^. Interestingly enough, 
the sum of these two values is very close to the cold molec- 
ular gas surface density we measured in the upper parts of 
the shock region. 

For the extended emission, the measurement of the cor- 
responding neutral and molecular gas masses cannot be re- 
alistically performed because the extended emission cover 
irregular parts of the large HI distribution in SQ and 
there are no CO observations covering the whole group 
with enough sensitivity to detect a plausible extended 
molecular hydrogen distribution. The measurement of the 
corresponding X-ray emitting gas mass is instead rather 
straightforward since it can be derived by the X-ray Halo 
luminosity, as for the shock region (see Appendix . Di- 
viding the total X-ray gas mass by the projected X-ray 
halo emission area (a circle of radius w 40 kpc), we deter- 
mined the hot gas mass surface density shown in col. 11 
of Table H 

7. DISCUSSION 

If one were to view SQ at a greater distance such that the 
group would appear as a point-like source to Spitzer (i.e. 
at a redshift of > 0.5 — 2, in the main star-forming epoch 
of the Universe, and also the epoch when galaxy groups 
were first forming) , one would not regard this as a partic- 
ularly unusual infrared source. On the basis of the shape 
of the total dust emission SED, plotted as the black curve 
in Fig. 1191 the only noteworthy points would be the quite 
warm MIR colours and moderately high FIR luminosities 
which would most likely lead our hypothetical observer to 
conclude that this source had an AGN, possibly combined 
with a mild starburst. In this he would be at least in part 
correct, as illustrated in Fig. jT^jby the curves represent- 
ing the constituent emission components from the AGN 
galaxy (red), from star formation regions (blue) and the 
emission from X-ray emitting regions (green), where one 
can immediately see the high relative contributions of the 
AGN galaxy to the 8/xm and 24/Lim emission, respectively 



~ 50% and - 70%. In the FIR, ahhough the AGN is 
still the most luminous individual source, there is emis- 
sion at a comparable level from the combination of the 
distributed star-formation regions in SQ and the infrared 
counterparts of the X-ray emitting shock and halo struc- 
tures; the combined FIR/submm SED resembles that of 
star-forming galaxies, with an amplitude similar to that 
of the local starburst galaxy M82. 

However our analysis of the spatially resolved structures 
has shown that the characteristics of infrared emission in 
SQ are the very opposite of a nuclear starburst, with the 
star-formation activity enhanced in regions far away from 
the main bodies of the galaxies. In the following Sect. 
17.1! and 17.2! we discuss to what extent the nature of this 
distributed star formation in SQ may differ from the star 
formation in the disks of individual galaxies, in terms of 
sources of gas and the star-formation efficiency for the 
group as a whole, and consider the related issue of the ex- 
tent to which collisional heating of grains in the IGM of 
SQ may be cooling the IGM and thus contributing to the 
fuelling of the star-formation. 

7.1. Star formation in SQ 

The results of the SED fitting to the constituent compo- 
nents of SQ indicate that star formation activity and pho- 
tons from old stars are the major agents powering the ob- 
served global dust emission from the group, supplemented 
by photons produced by the accretion flows in the AGN 
and possibly, in the case of the X-ray emitting regions, by 
collisional heating. In our quantitative discussion of star 
formation activity in SQ we will adopt an initial working 
hypothesis that the collisional heating mechanism is mi- 
nor compared with photon heating in the X-ray sources. 
We will scrutinized this hypothesis in detail in Sect. 17. 2i 
where we discuss physical constraints of the fraction of the 
infrared emission that can be collisionally powered. 

Under this working hypothesis we can gather together 
the information from Table H| to obtain a total global star 
formation rate of 7.5 MQ/yr for SQ in its entirety. This 
global star formation rate does not seem particularly dis- 
crepant from that typically found for most of the galaxies 
of similar mass in the local Universe that are clearly not 
so strongly interacting as SQ galaxies. Using the empiri- 
cal relation between galaxy stellar mass and SF R for field 
galaxies shown in Fig. 17 of iBrinchmann et al.l f2004ll we 
can estimate the typical SFR of galaxies having the same 
stellar mass of SQ galaxies, of order of w 10^^ Mq (see 
table |S|). For this value of the stellar mass, the mode of 
the distribution is at SFR k, 1 MQ/yr. Therefore, since 
we estimated the SFR in a field containing four galaxies, 
the expected SFR would be w 4MQ/yr, comparable to 
the measured value. Thus it seems that the star forma- 
tion efficiency of SQ in relation to field galaxies is largely 
independent of whether the gas is inside or outside the 
main stellar disk. In fact, of the total global star forma- 
tion rate of 7.5 M0/yr for SQ, 2.2 and 5.3 M0/yr can 
respectively be ascribed to the SED components for star 
formation regions and X-ray sources in Fig. 1191 This is 
a very remarkable result indicating, as it does, that the 
bulk of star formation activity in SQ is apparently asso- 
ciated with X-ray emitting structures, occurring far away 
from the galaxy centers, either at the peripheries of the 
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galaxies or in the intergalactic medium. Furthermore, the 
total extragalactic SFR is well in excess of the SFRs of 
the previously studied individual examples of extragalac- 
tic compact star formation regions, SQ A and SQ B. 

To quantify the efficiency of star formation in the var- 
ious components of SQ in relation to that of the disks of 
individual galaxies, we plotted in Fig. [^Hlthe SFR per unit 
area for dust emission sources in SQ as a function of the 
gas mass surface density. On the same diagram, as ref- 
erence , we plotted the relations found by iKennicutt et al.l 
(|2007l) (K07) and iBigiel et al.l EOO^ (BOS) for star for- 
mation regions inside nearby spiral galaxies. This figure 
shows a wide range of star formation efficiencies, which 
we discuss below for each of the structural components, 
with special emphasis on the extended component of star 
formation which dominates the global SF activity in SQ. 

7.1.1. Discrete Star Formation Regions 

Fig. EOl shows that SQ A and SQ B, have star formation 
rates very similar to those observed for galaxian regions 
with the same gas column density, despite their being lo- 
cated well outside the galaxies of the group . This result 
is con sistent with the conclusions reached bv lBraine et al.l 
(|2001|). studying a sa mple of tidal dwarf galaxies, and by 
iBoauien et al.l ( 20091) . which performed a multiwavelength 
analysis of star formation regions in collisional debris. The 
situation is however different for the other two bright star 
formation regions, HII SE and HII SW, which present 
much higher SFRs than those predicted by the plotted 
relations. This might in principle be due to a more effi- 
cient mode of star formation happening in those regions. 
It would be of interest to acquire more complete infor- 
mation on the gas content of these sources to further in- 
vestigate this conjecture. An alternative explanation for 
the high SFR found in HII SE and HII SW is that there 
is an additional component of UV emission produced by 
the radiative cooling of shocked gas. These UV photons, 
unrelated to star formation phenomena, can in principle 
contaminate the UV flux measurement but also power part 
of the observed dust emission, thus leading to an overesti- 
mation of the SFR. Optical spectra of these regions show 
evidences of line shock excitation (P.-A. Due and collabo- 
rators, private communication), therefore one cannot rule 
out this possibility. 

7.1.2. Star Formation associated with the Shock 

The average SFR in the shock region (3.1 x 
10~^ M0/yr/kpc^) seems to be well in agreement with 
the empirical SFR-gas surface density relation in Fig. [201 
Thus, at first sight, the shock does not seem to have had 
much effect on the star formation activity in the stripped 
interstellar gas, neither triggering enhanced star forma- 
tion through shock-compression of dense clouds of gas, 
nor suppressing star formation through heating and dis- 
persion of the clouds. The observational situation is how- 
ever complex in that the star formation observed towards 
this region could be happening inside the extended fea- 
tures connected with the incoming galaxy NGC 7318b, 
seen through, but not obviously physically co-existent with 
the shock region. This is supported by the local morphol- 
ogy of the UV and MIR emission that seems to follow 
the optical shape of the intruder, instead of showing a 



linear north-south ridge resembling the shocked gas emis- 
sion morphology. These two scenarios could in principle be 
combined since the ISM of the intruder galaxy has presum- 
ably been shocked as well and this could have triggered the 
observed star formation out of the gas associated with the 
intruder galaxy. In addition, it is conceivable that some 
fraction of the UV luminosity, albeit probably a minority, 
could be due to gas cooling rather than from massive stars. 

It is also of interest to compare our constraints on the 
SFR in the shock region with the luminosity of the radio 
synchrotron emission from the shocked gas. Traditionally 
radio synchrotron measurements are compared to infrared 
emission measurements in the context of the radio-FIR 
correlation for indivi dual galaxies , for w hich we use here 
the relation given bv lPierini et al.l ()2003t ): 

log Li.4 GHz = 1.10 X log Lfir - 18.53 (7) 
where L1.4 qhz is the luminosity at v = 1.4 GHz in units 
of W Hz^^ and LpiR is the total FI R luminosity i n units 
of W. We used the relation from iPierini et al.l (|2003() 
since it was derived from data covering cold dust emis- 
sion longwards of > 100/im, where, as in the case for 
the SQ shock region, most of the power is radiated. Us- 
ing the value we measured for the shock region FIR lu- 
minosity, Lfir = 1.5 x lO'^^W, we obtain for the pre- 
dicted radio luminosity L1.4 ghz = 1-8 x 10^^ W Hz. 
This value is 20 times smaller than the radio luminos- 
ity derived from the radio measurement in Xu et al. 03: 
Li,4 GHz = 3.7 X 10^2 W Hz. The high radio/infrared lu- 
minosity ratio (which would be higher still with respect 
to the photon-powered component if collisional heating 
was important) points to the existence of an additional 
source of relativistic particles in the SQ shock, accelerate d 
at the shock itself (see e.g. iBlandford fc EichleJ 119871 ). 
which dominates the population of particles accelerated in 
sources more directly linked to star formation regions such 
as supernova remnants. This confirms and strengthens 
the preliminary results of lXu et"al] ()2003[ ). and, bearing in 
mind that the total radio emission from SQ is dominated 
by the emission from the shock region, suggests that cau- 
tion is needed in using radio synchrotron measurements 
to infer star formation rates in groups involving strong 
dynamical interactions of galaxies with the IGM. 

7.1.3. Star Formation associated with the Extended 
Emission 

Our measurements have shown that both the obscured 
and visible components of the star formation in SQ are dis- 
tributed in a widespread pattern, loosely coincident both 
with the overall dimensions of the group as well as with 
the extended "halo" of X-ray emission. This current mor- 
phology of SQ star formation can readily be accounted for 
in terms of the several galaxy - galaxy or galaxy - IGM 
interactions which have occurred in the group. In this 
standard scenario, the gas that is currently converted into 
stars is interstellar in origin (as we will discuss later in 
this section). However, the data are also consistent with a 
rather different scenario, in which the reservoir of gas out 
of which stars are being formed is the hot gas phase of the 
intragroup medium. 

Specifically we consider a scenario in which the hot in- 
tergalactic medium is cooling and condensing into clouds 
which are the sites of the star formation providing the UV 



12 



Natale et al. 



powered component of the extended infrared emission. If 
this is occurring in a steady state, in which the rate at 
which the removal of hot gas by coohng is balanced by ac- 
cretion of further primordial gas (that is, we assume that 
variation in the accretion rate have longer time scales than 
the cooling time scale) , we can write 



Thot 

where Mhot and rhot respectively denote the mass and 
cooling timescale of the X-ray emitting medium, and the 
inequality denotes the fraction of the cooling gas that 
ultimately condenses into stars. If the observed X-ray 
emission is the main component of luminosity of the hot 
medium, and taking observational values for the X-ray 
emitting plasma from Trinchieri et al. (2005) of A/hot — 
1.0 X lOi°M0 and Thot » 0.3 x W^^yi (see Sect. 7.2.2) we 
obtain S'Fi?hot < SMQyr^^, comparable to the observed 
value of 5.3MQyr~^ for star formation seen towards the X- 
ray-emitting halo. In reality, however, one would expect a 
large fraction of the cold gas resulting from the cooling of 
the gas to be recycled into the hot medium due to the feed- 
back of mechanical energy from the newly-formed stars, 
in which case SFRi^ot would be substantially less the up- 
per limit from Eqn|Sl Thus, if X-ray emission driven by 
gas-gas collisions is the dominant cooling path for the X- 
ray halo, one would conclude that the bulk of the star 
formation would not be directly fuelled out of the IGM. 
On the other hand, as shown in Sect. 7.2.2, our Spitzer 
FIR data, coupled with astrophysical constraints on the 
injection rate of grains into the hot medium do not at 
present rule out the possibility that the dominant cooling 
mechanism of the halo is FIR emission driven by gas-grain 
collisions. At present, therefore, we cannot rule out star 
formation out of a primordial IGM purely on considera- 
tions of gas fuelling. 

A more powerful constraint would be to consider the ef- 
ficiency at which such a mode would have to operate at. 
Unfortunately, the total mass of molecular cold gas, A/coid 
on a global scale in the halo is unknown, preventing a di- 
rect empirical measurement of star formation efficiency of 
the halo on the KS diagram. One can however estimate 
A^coid under our simple steady state scenario by writing 

Thot Tcold 

where Tcom is the typical time scale spent by the gas in 
the cold molecular phase before being converted into stars, 
which is the timescale for the collapse of a molecular clouds 
to form stars, multiplied by the mean number of times the 
cold gas is recycled into the hot medium through mechan- 
ical feedback before condensing into a star. Even allowing 
for several cooling cycles of the hot gas before condensing 
into stars, the very short collapse timescale of a few million 
years for molecular clouds - some three orders of magni- 
tude shorter than Thot - makes it likely that Thot > > TcoW , 
which in turn implies Mhot >> Mcoid- If this is the case, 
the cold gas amount, related to the extended emission, that 
we should add to Sgas would be negligible, and the position 
of the point for the SQ halo on the KS diagram would be 
well to the left of the relation for star-forming galaxies, as 
shown in Fig. 1201 This would require a more efficient star- 
formation process in the IGM than is typically observed 



to occur in the disks of star-forming field galaxies, which 
would seem to be counter-intuitive to the naive expecta- 
tion that denser, colder reservoirs of gas closer to the min- 
imum of the gravitational potential wells associated with 
individual galaxies should form stars more easily. On the 
other hand a definition of efficiency in terms of gas surface 
density over scales much larger than the star-formation 
regions may have limited predictive power in this context. 
Furthermore, any star formation in a cooling hot IGM 
would not be expected to occurring in rotationally sup- 
ported systems similar to the galaxies used to define the 
KS relation, so may have a different relation to the local 
gas density. 

The second, more conventional scenario to explain the 
widespread star formation in the extended halo compo- 
nent of SQ, is to invoke interstellar gas as the source of 
cold gas fuelling most of the star formation. In this sce- 
nario tidal interactions between galaxies or hydrodynam- 
ical interactions of the ISM with the IGM remove cold 
galaxy material that, for SQ, produces a similar SFR as 
would have been the case for a similar amount of gas inside 
field galaxies. For groups in general the statistical relation 
of neutral gas observed within and outside the member 
galaxies to the X-ray emission characteristics of the IGM 
indeed indicates that at least some of the cold gas in the 
IGM medium originates in the galaxies, and that at least 
some of the X-ray emitting gas also has an interstellar ori- 
gin ( eg. IVerdes-Montenegro et al.lF2001l : iRasmussen et al.l 
|2008[) . More direct evidence of the removal of interstellar 
gas and dust through galaxy interactions in high density 
environments is provided by observations of interacting 
galaxies in clusters, for example through Herschel imag- 
ing of extraplanar dust emission associated with stripped 
atomic and molecular gas around the Virgo cluster galaxy 
NGC4438 (Cortese et al, 2010). 

In SQ, there is direct morphological evidence for burst- 
ing sources associated with tidally removed interstellar gas 
in the form of SQ A and SQ B, and a hint that other less 
prominent components of the FIR emission in the IGM 
may also be associated with tidal features, such as the 
enhancement of Ha and MIR/FIR emission seen towards 
the bridge linking the intruder galaxy with the AGN host 
galaxy. In this picture the extended emission is a conglom- 
eration of discrete sources similar in nature to SQ A and 
SQ B but much more numerous and of lower power. The 
fact that SQ A and SQ B are likely to have a fundamen- 
tally interstellar origin despite their present location out- 
side the galaxies is further supported by the rather high 
value of the dust to gas ratio we have found for these 
sources, Zd — 0.002, which is a signature of a high gas 
metallicity. The one puzzling aspect of the extension of 
this scenario to the global star formation in the IGM of SQ 
is the lack of prominent star formation associated with the 
bulk of the HI in the IGM located to the So uth and East of 
the group. However, as shown by Fig. 5 of IWilliams et al.l 
(200^, the gas column density in these clouds (generally 
less than w 3 x 10^° atoms cm~^) is rather low compared 
to the values found for the brightest star formation sites 
in SQ. Therefore it is plausible that star formation is in- 
hibited in these clouds because of the low gas density. 

Finally, we remark that if stripped neutral and molecu- 
lar interstellar gas was the only source of cold gas powering 
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star formation, one would expect a rapid quenching of star 
formation activity in SQ in the future, whereas if the star 
formation were fuelled by a cooling primordial IGM, the 
star formation activity would simply follow the accretion 
rate onto the group of the primordial IGM. Thus, mea- 
surement of the UV and MIR/FIR luminosity functions of 
well defined statistical samples of galaxy groups, now be- 
ing defined through deep optical spectroscopical surveys 
such as the Galax y And Mass Assembly survey (GAMA; 
[Driver et al1l2009D . may offer a way of quantifying the rel- 
ative importance of these modes of star formation on the 
dynamical halo mass of groups in the local Universe. 

7.2. The nature of the X-ray correlated FIR emission in 
Stephan's Quintet 

Previous searches for infrared emission counterparts of 
the hot X-ray emitting components of the intergalactic 
medium of nearby objects have almost exclusively targeted 
the intraclustcr medium (IGM) of rich galaxy clusters, ei- 
ther using a stacking analysis for X-ray or optical selected 
clusters from the IRA S all sky survey (iGiard et al.ll2008t 
iRoncarein eralll2010 ) or in detailed imaging observations 
of the Goma Cluster and other clusters with ISO or Spitzer 
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iBai et al.l l2007riKitavama et al.H 2009'). Ah these studies 
have either yielded upper limits or marginal apparent de- 
tections of the IGM at far-IR brightness levels far fainter 
than those we have measured towards the X-ray emitting 
IGM of SQ with Spitzer. Furthermore, even the apparent 
detections of the IGM were susceptible to confusion with 
foreground cirrus or the back ground galaxy popula tion due 
to the accidental similarity (jPopescu et al.l [2000() in far- 
IR colour of the collisionally- and photon-heated emis- 
sions, so realistically must be treated as upper limits to 
any c ollisionally heated emission from the IGM (|Bai et al.l 
l2007f ). Although intracluster dust has been unambigu- 
ously detected at optical wavelengths, through statistical 
studies of the reddening of bac kground sources through 
large numbers of clusters (Ghelo uche et al.|l2007[ ) this has 
been at a low abundance of about 10"'' in the dust-to-gas 
mass ratio, consistent both with the non-detections of the 
IGM in infrared emis sion and with specific predictions by 
iPopescu et al.l ()2000t ) of the dust content of typical IGMs. 

In the present work we have been able to sidestep the 
confusion problems afflicting the previous attempts to de- 
tect infrared emission from the IGM in clusters only by 
virtue of the much higher far-IR brightness levels mea- 
sured towards the X-ray emitting structures in SQ. These 
high brightness levels, together with the correspondence of 
the infrared emission morphologies with those seen in the 
UV, and the evidence from our fits to the infrared emis- 
sion SEDs described in Sect. 15.31 however suggest that a 
major part of the dust emission is photon heated. This, in 
turn, raises an apparent paradox that the closest spatial 
correlation between dust emission and gas column density 
is not, as might be expected for photon heating, with the 
cool HI and molecular gas component, but rather with the 
hot X-ray emitting component. In our discussion on star 
formation in SQ we identified a possible way out of this 
paradox by postulating that the extragalactic star forma- 
tion in SQ is fuelled by gas from a hot IGM, whose cooling 
is enhanced by grains injected into the hot IGM. Here, we 



use the relative levels of detected emission at infrared, X- 
ray and UV/optical wavelengths to constrain the possible 
sources of grains and the resulting collisionally driven cool- 
ing of the hot plasmas, considering separately the two main 
X-ray emitting structures corresponding to the shock and 
the halo regions. 

7.2.1. The shock region 

The main constraint on the coUisional heating of dust 
in the X-ray emitting plasma downstream of the shock is 
given by the high dust-to-gas mass ratio that would be 
required to reproduce the FIR measurements in the case 
that grain heating by a diffuse radiation field can be ne- 
glected. For pure collisionally heated emission from grains 
in a homogeneous medium, a dust mass Afd = 2 X 10^ M0 
(see Sect. 15. 3p would be required, which, taken together 
with the measured X-ray emitting gas mass, Mx « 10^ Mq 
(see Appendix [C]), would imply a dust-to-gas ratio for the 
hot plasma of Zd = 0.02. This value is high in comparison 
with expectations based on a balance between injection of 
grains and their removal through sputtering. Assuming 
that the shock is propagating into a stripped interstellar 
medium or through the ISM of the intruder galaxy, char- 
acterized by a dust-to-gas ratio comparable to the Milky 
way value, we can derive the amount of dust mass per unit 
time A/d that is injected downstream of the shock: 

Md = VskAp^pZd,np (10) 

wher e 14h = 600 km/s is the shock velocity (jGuillard et al.l 
I2009D . A = 330 kpc^ is the projected area of the shock 
front (corresponding to 20 x SOarcs^ on the sky), p^p « 
l.lriHTOH = 4.5 X 10~^^ g/cm^ is the upstream gas den- 
sity (riH = 0.01/4. cm""^) and Z^^up = 0.01 is the assumed 
dust to gas ratio for the gas upstream of the shock front. 
The resulting dust injection rate is A^d = 0.14 M0/yr. To 
estimate grain destruction through sputterin g in collisions 
with heavy particles we use the formula of iDraine et al.l 
()1979() for the dust sputtering time scale: 



' sp 



2 X 10^ — yr 



(11) 



where a is the dust size in pm and tih is the proton num- 
ber density in cm~^. Assuming tt-h ~ 0.01 cm"'^ and 
a = 0.1 /im, we get Tgp = 2 x 10^ yr. Explicit calculation 
of the cooling of the shock gas due to gas-gas and gas- 
grain collisions, taking into account the evolution of the 
grain size distribution (under the assumption that there 
are no further sources of grains) , show that the gas cool- 
ing time scale is icooi ~ 7 x lO'' yr, more than three times 
larger than r^p (Natale G.,2010, PhD Thesis). Under these 
circumstances a conservative estimate for the dust to gas 
ratio downstream of the shock can be derived by dividing 
the total mass of dust predicted to be at any time down- 
stream of the shock to the measured hot gas mass, that 
is: 

MdTsp 



■^d.down — 



(12) 



Applying this formula, the expected value for the X-ray 
plasma dust to gas ratio is ^d,down ~ 3 x 10"'^, a fac- 
tor w 7 smaller than the value inferred from the data, for 
the case where the entire FIR dust emission is powered 
by coUisional heating. This is only a crude estimate for 
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the dust abundance, which might take higher values if the 
hne of sight depth of the shock was assumed to be larger 
than the projected width seen on the sky, but which on 
the other hand is an overestimate compared to predictions 
taking into account the true evolution of grain size down- 
stream of the shock. In any case, it is evident that the 
measured dust to gas ratio is much larger than the pre- 
dicted 2'd.down, and a further source of dust grains would 
be required if collisional heating were to be a significant 
driver of the FIR emission. 

The most obvious such source would be the reservoir 
of grains in the cold component of the gas in the shock 
region, which could be potentially released into the hot 
gas by ablation if the clouds were in relative motion to 
the hot gas. As already mentioned the presence of cold 
gas is implied by the detection of 8/im emission and has 
been directly demonstrated by the de tection of molecula r 
hydrogen in th e MIR rotational lines (jCluver et al.ll2010D . 
iGuillard et al.l ()20Q9[ ) modelled the H2 emitting clouds as 
cooling pre-existing clouds in the upstream medium, a sce- 
nario which would predict the clouds to be in relative mo- 
tion to the hot gas downstream of the shock. Another pos- 
sibility which might potentially account for a FIR counter- 
part to the shock is to invoke heating of grains upstream 
of the shock by the UV radiation from the cooling clouds 
downstream of the shock. A quantitative treatment of all 
these effects would however require a self consistent model 
for the passage of the shock through a two phase medium 
which tracked the exchange of photons, gas and dust be- 
tween the phases, which is beyond the scope of this paper. 

It's appropriate to emphasize here that even in the case 
that the collisional heating is relatively unimportant in 
bolometric terms, the effect on the gas cooling time scale 
is not negligible. Detailed calculations for a steady state 
ID homogeneous model, considering dust sputtering and 
dust cooling, show that the cooling time scale is shorter by 
a factor « 2 - 3 (Natale G.,2010,PhD thesis) for shock ve- 
locities and gas densities similar to those characteristic of 
the shock in SQ. From this perspective it would be impor- 
tant to have more direct empirical constraints on the gas 
cooling mechanism using the improved sensitivity, angular 
resolution and wavelength coverage of the Herschel Space 
Observatory. Using the simple approach described above, 
we estimate the surface brightness of the radiation emitted 
by coUisionally heated dust to be equal to 1.4 MJy/sr at 
150/xm, the wavelength where the flux density is predicted 
to peak according to the collisional heated dust model we 
used for the shock region SED fit (see Fig. [TS]) . 

7.2.2. The extended emission 

As for the shock, the detection of substantial 8/im emis- 
sion towards the extended X-ray emitting halo indicates 
the presence of a component of dust in cold gas phases, 
shielded from collisions with electrons and ions in the hot 
plasma, which can only be heated by photons. Rather 
interestingly, the curves of growth we constructed after 
removing or masking the galaxies and the compact star 
formation regions from the infrared and FUV maps (Fig. 
[TU|) show a very similar profile at all wavelengths. This is 
a strong indication that FUV sources are associated with 
the FIR extended emission and, therefore, star formation 
is the main mechanism powering the extended emission. 



The total luminosity of the extended dust emission is 
4 X lO'*'^ erg/s, a factor ss 130 higher than the X-ray Halo 
luminosity (Trinchieri et al.2005, see also Appendix \C\\ . 
This large difference between the dust emission and X-ray 
luminosities implies that collisional heating of dust grains 
could still be important for the cooling of the hot halo gas 
even if it only powers a small fraction of the dust luminos- 
ity. We define the fraction i^coU of dust emission that is 
coUisionally powered, by the relation: 

-^ir" = -FcouilR (13) 
where if^ ' is the dust luminosity that is powered by col- 
lisional heating. 

A first physical constraint on i^coii that we can apply is 
the requirement that the gas cooling time scale TcooI is not 
larger than the dynamical hot gas time scale, given by the 
sound crossing time Tcross- Including dust cooling, TcooI is 
equal to: 

3 kT ( Mx \ 
Tcooi = 77 1 ; — ) (14) 



2 ^mn V -^x + FcowLyr. J 
where T — Q x 10^ K is the gas temperature, Afx = 
10^° Mq is the X-ray emitting gas mass, and Lx and Lir 
are the X-ray and infrared luminosities. Requiring that 
T'cooi > Tcross = 2i?haio/cs, whcrc i? = 35 kpc is the halo 
radius and Cg ~ 300 km/s is the sound speed, one obtains 
^^coii < 0.05. 

Another constraint can be derived from an estimation of 
the dust injection rate into the hot gas. For equilibrium 
between dust injection and sputtering, one would have: 

Md^i^coii— = 2.4Feon ^ (15) 
Tsp yr 

where is the dust injection rate, Tgp 2 x 10^ yr is 
the sputtering time scale (calculated using formula lTll with 
UH = O.OOlcm-3 and a = O.l^im), and Md = 4.7 x 10^ Mq 
is the inferred dust mass producing the extended emission. 

The first dust injection sources we consider are the halo 
stars. As described in Appendix |E1 we have derived an up- 
per limit to the dust injection rate from halo st ars based on 
the o bserved R-band halo surface brightness (| Moles et al.l 
Il998f ) and the theoretical predictions for Stardust injec- 
tion by Zhukovska et al.OS: Md = O.O75M0/yr. Substi- 
tuting this value in Eqn. [T5l we get i^coU < 0.03. A 
similar procedure could be followed considering the dust 
injection from type II supernovae given the inferred SFR 
associated with the extended emission. Empirical stud- 
ies of Galactic supernova remnants have shown the yield 
of condensate pe r supernova event to be of the order o f 
0.01 to 0.1 M0 (|Green et al.ll2004 iFischera et al.ll2002[ ). 
Relating the frequency of core collapse supernova events 
to the SFR in the extended component of SQ, we esti- 
mated Md = 7 X lO"'^ Mg/yr assuming 0.01 Mq of dust 
produced in each supernova. Therefore supernovae dust 
injection cannot be considered important in replenishing 
dust in the hot med ium. Finally w e considered the mech- 
anism proposed by iPopescu et al.l ({2000.) where dust in 
IGM external to a cluster can be introduced into the hot 
intracluster medium in the supersonic accretion flow of 
barions onto the cluster. In the case of SQ the baryonic 
accretion rate can be estimated to be of the order of the 
current IGM seen in X-ray divided by a Hubble time which 
is ?a 2 Mq /yr. Even if we consider that there is maybe a 
cooler comp onent of the accr eting IGM not visible in the 
X-ray (|Mulchaev et al.l [20001) . it does not seem plausible 
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to achieve dust accretion rate for the particular case of SQ 
much greater than 0.01 M0/yr for reasonable values of the 
dust to gas ratio in the accreting material. 

Given these estimates, we can conclude that if there is 
a contribution of coUisional heated dust to the total ob- 
served dust emission, this should only be of the order of 
a few percent and maintained by direct injection from the 
in situ stellar population. This however is sufficient to 
maintain a coUisionally powered FIR luminosity of up to 
PS 10^^ erg/s, which is several times the X-ray luminosity. 

Such a level of FIR cooling could help to provide a nat- 
ural explanation for the rather low ratio of ~ 0.02 for 
the ratio of the mass of baryons in the hot X-ray emit- 
ting halo of Mx = 2 X 10^° Mg to the dynamical 
mass of ^ 10^^ M^P^. For the concordance CDM cos- 
mology, this ratio should be ~ 0.09 in the ab sence of sig- 
nifica nt dissipative effects (see e.g. Fig. 9 iFrenk et al.l 
Il999l ). Recalling that the total baryonic mass currently 
locked in stars in SQ of ~ 10^^ Mq is comparable to the 
expected total mass of baryons, it is apparent that this 
scenario is also quantitatively consistent with the efficient 
condensation of the IGM accreting onto the group into 
stars, as postulated in Sect. 7.1. Some empirical indica- 
tio n of efficient star formation in groups has been found 
bv lTran et al] (|2009D . who showed that the group environ- 
ment is more conducive to star formation than either the 
cluster or field environment. In general, enhanced cool- 
ing due to a FIR coUisionally powered luminosity of a few 
times Lx would, since it promotes the condensation of hot 
X-ray emitting gas into stars, also provide an alternative 
explanation for the steeper Lx — Tx relation observed for 
groups with lower virial velocities, an effect that is other- 
wise attributed to thermal or kinetic feedback from AGNs 
(e.g. ICavaliere et"aIll2008D . 

In summary, the detection and recognition of the puta- 
tive coUisionally heated component of FIR emission from 
SQ will be a crucial measurement with wide ranging im- 
plications for the nature of star-formation in the group, 
the division of Baryons between hot gas, cold gas and 
stars, and the thermodynamic properties of the IGM in 
SQ. Such a detection will however be very challenging, re- 
quiring a combination of excellent surface brightness sen- 
sitivity and angular resolution at submm wavelengths, as 
illustrated by the following estimates: assuming the upper 
limit Fcou = 0.03, we predict the SED of the coUisionally 
powered component of the extended emission to peak at a 
level of 0.15 MJy/sr at 190/^m. 

7.3. The AGN galaxy NGC 7319 dust emission 

The Seyfert2 galaxy NGC 7319 is the most powerful 
source of IR emission in SQ and is the only galaxy where 
we found enhanced infrared emission in the central regions. 
The total IR luminosity is 9.6 x 10'*'^ erg/s. Of this, we 
deduced from the SED fitting using a combination of the 
AGN torus template and the diffuse dust emission compo- 
nent (see Sect. 15. 2p that 4.1 x 10^'^ erg/s are emitted by 
dust in the torus and 5.5 x 10**^ erg/s by diffuse dust. As 
one can see from the SED best fit, shown in Fig. [131 the 

This value is derived from the X-ray temperature that is consistent with the velocity dispersion of the group, a = 340km/s, as calculate d by 
lO 'Sullivan et af] ^00§). However including the more recent recession velocity measurement for NGC 7320c, quoted in Sect. 4.1 of Sulontic et al.l 
Il200in . the group velocity dispersion would he (t = 80km/s indicating a DMH mass of order of IO-'^-'^Mq. Since the number of group galaxies is 
very low for a good estimate of the velocity dispersion, we consider the X-ray temperature based estimate for the DMH mass more reliable in 
this case. 



24/im MIR emission is dominated by the torus emission 
which also contributes « 25% of the 8/Ltm emission. This 
is consistent with the observed morphology at 8 and 24/im 
emission, which is a combination of a central source plus 
a disk component in the case of the 8/xm map and a dom- 
inant central source in the case of the 24/im map (see Fig. 
[1]). According to the SED fit, the FIR 70 and 160/im emis- 
sion is due to diffuse dust emission. This is also consistent 
with the brightness distributions on the Spitzer FIR maps 
which show the emission to be partially resolved, with ex- 
tents inferred from the FIR map fitting technique of 5 
and 12 kpc at 70 and 160/im respectively. This extended 
FIR emission is centred at the position of the AGN, po- 
tentially indicating the presence of star formation regions 
distributed on kpc scales in the su rrounding areas. How- 
ever, as one can see from Fig. 1 of iGao et al.l (l200l . the 
CO emission is predominantly from the northern regions 
of this galaxy. Therefore the locations of the bulk of cold 
gas, where the star formation sites should be, and the cold 
dust emission are not coincident suggesting that star for- 
mation may not be the main phenomena powering the FIR 
dust emission. 

Overall, the diffuse dust emi ssion in the inner disk , to- 
gether with the absence of HI ([Williams et al.|[2002[ ) and 
the displacement of the molecular hydrogen emission sug- 
gest that galaxy-galaxy interactions have removed gas 
from the central regions on kpc scales, as has plainly also 
happened in the other late-type SQ galaxies. In the case 
of NGC 7319, the interactions have apparently also in- 
duced gas flows that have triggered AGN activity. This 
raises the question whether the AGN itself might be the 
source of photons powering the surrounding diffuse dust 
emission in this galaxy. According to the standard AGN 
unification theory, the AGN torus of a Seyfert 2 galaxy is 
viewed edge on, which, in the case of NGC 7319, would 
imply that the plane of the torus is almost perpendicular 
to the galaxy disk, since NGC 7319 is seen basically face 
on (see SDSS g-band map in Fig. [T]). In this configura- 
tion, it would indeed be plausible that accretion powered 
radiation, escaping into the polar direction from the AGN 
torus structure could propagate into the galaxy disk and 
contribute to the diffuse radiation field in the central parts 
of the galaxy. An especially efficient coupling between the 
dust and radiation would be expected if the polar axis of 
the AGN torus were to be parallel to the diffuse dusty disk 
of the host galaxy. 

One indirect indication that the radiation from an AGN 
with this orientation may be responsible for at least part 
of the diffuse FIR emission in NGC 7319 comes from the 
rather low observed UV/FIR luminosity ratio, which is 
qualitatively consistent with the expected strong atten- 
uation of the direct UV light from an AGN to which the 
line-of-sight passes through the torus. From the measured 
FUV luminosity, Lfuv = 4 x 10**^ erg/s, and from the 60 
and 100/im luminosities extrapolated from the SED best 



fit, L 



60/j,m 



2 X 10^3 gj.g/s and L 



lOOAim 



3.7 X lO'*^ erg/s, 

we have found LFUv/^eoAim = 0.2 and Lpuv/iiooMm = 
0.1. These values are several factors lower than the av- 
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erage values found by iPopescu et alj ()2002[ ) for galax- 
ies of similar morphologic al type (L\w / L ^num = 1 and 
Ltjy / LiQOfim — 0.6) and bv lXu et"aLl ()2006[ ) for galaxies of 
similar iFUV+^^eo/im total luminosity (i'uv/-t'60Aim = 0.5). 

A scenario in which the diffuse dust emission in NGC 
7319 is powered by the AGN also seems to be broadly 
consistent with simple constraints on the UV luminosity 
of the AGN derived from the observed hard X-ray lu- 
minosity. Specifically, for L2-10 keV = 2.8 x 10^^ erg/s 
(|Trinchieri et al.ll2005l) and adopting the bolometric cor- 
rections for AGN gal axies with similar L2-10 keV by 
IVasudevan et all (|2010D . 

'«2-io keV ~ 10 — 30, one pre- 
dicts a total UV luminosity from the AGN of Xagn = 
K2-10 kcvi2-io kov = 2.8 - 8.4 X 10^^ erg/s. The upper 
end of this range is comparable to the total MIR/FIR lu- 
minosity of the AGN and AGN host galaxy. If the AGN 
were indeed powering the diffuse dust emission in the in- 
ner disk one would expect that future FIR observations 
of angular resolution sufficient to image the diffuse dust 
emission could be used to determine some basic geometri- 
cal properties of the AGN. Specifically, one would expect a 
cone structure for the component of diffuse dust heated by 
the AGN, with a strong radial colour gradient, and with a 
symmetry axis aligned with the polar axis of the torus. 

Finally, we draw attention to a further FIR-emitting 
feature potentially connected to the AGN, whose origin 
has still to be clarified. This is the bridge which appar- 
ently connects the AGN to t he center of the shoc k re- 
gion, ap pcarmg m the X-ray ()Trinchieri et al.ll2005D. Ha 
(IXu et al. 1999) and warm H2 line emission ( Cluver et al.l 
l201Clf ). The presence of a cold dust emission from this fea- 
ture is hinted by a slight asymmetry of the FIR emission on 
the 160/Ltm residual map (see Fig. 0]). High resolution FIR 
observations would also be useful to elucidate the corre- 
spondence between dust emission and the gas distribution 
in this bridge feature. 

8. SUMMARY 

In this work we have performed an extensive analysis of 
the infrared emission from SQ, based on the Spitzer MIR 
and FIR maps and including multiwavelength data from 
optical/UV, Radio and X-ray. The main aim has been un- 
derstanding the nature of the FIR emission seen on the 
Spitzer 70 and 160/im maps and its implications for star 
formation and hot gas cooling in this group of galaxies. 
The technical part of our analysis can be summarized in 
four points: 

1) We have identified the FIR sources seen on the Spitzer 
maps and performed their photometry on both the MIR 
and FIR maps. For the FIR photometry we devised a new 
source fitting technique that we used to model the source 
emission on the maps (see Figs. [2]and|3l). We measured 
MIR fluxes using aperture photometry (see Figs. [S] and 
[7]). Photometry results are shown in Tabled] 

2) We have compared the residual FIR emission (defined 
as the FIR emission after the subtraction of galaxies and 
discrete star formation regions) with HI distribution, X- 
ray emission and UV emission (see Fig. |4]). This compar- 
ison shows that the residual FIR emission, dominated by 
the emission from the shock region and a previously unde- 
tected extended component, is correlated with X-ray emis- 
sion and anticorrelated with the HI distribution. At the 



same time, the area covered by the FIR residual emission is 
populated by an extended distribution of UV sources asso- 
ciated with tidal features and intergalactic star formation 
regions. 

3) We have modelled the inferred infrared source SEDs 
using combinations of dust emission templates. These 
templates, described in Appendix [Bl are used to reproduce 
emission from PDR/HII regions, diffuse photon-powered 
dust emission, collisionally powered dust emission and in- 
frared emission from AGN torus. Among the immediate 
outcomes of the fitting procedure there are the dust mass 
required to produce the observed emission, the total dust 
infrared luminosity and the amount of UV powered dust 
emission (see Table [3]). 

4) Combining the amount of absorbed UV luminosity, 
obtained from the infrared SED fitting, with the observed 
UV luminosity, measured directly on the GALEX FUV 
maps, we determined the star formation rate for each dust 
emitting source. We also measured gas masses/column 
densities from the observed radio HI and CO lines and soft 
X-ray fluxes and, combining these measurements with the 
dust masses inferred from the infrared SED flts, we esti- 
mated the dust to gas ratios for several sources. Results 
are summarized in Table SI 

From this set of measurements and morphological com- 
parisons, together with simple predictions for the gas cool- 
ing time scale and the dust injection rate in hot plasma, 
we derived the following conclusions: 

1) The total star formation rate, within a circular area 
of radius 90" containing the four group galaxies in close 
proximity, is 7.5M0/yr. The star formation regions in the 
group are localized at the edges of the galaxies, on tidal 
features and in the intergalactic medium; no evidences for 
star formation have been found towards galaxy centers. 

2) The star formation efficiency for the star formation 
regions SQ-A and SQ-B, localized on tidal features, is close 
to that observed for star formation regions within nearby 
spiral galaxy disks (see Fig. [^U]) . Star formation in HII-SE 
and HII-SW seems to be much more efficient but the lack 
of complete gas measurements and the possible contribu- 
tion to UV radiation from radiative shocks does not allow 
a definite conclusion for these sources. 

3) The star formation efficiency in the shock region 
seems to be consistent with the star formation observed 
in nearby spiral galaxies. However, because of the elon- 
gated tidal feature of the intruder galaxy NGC 7318b seen 
through the shock region, it is not clear if the observed star 
formation in that area is associated with the shocked gas 
or not. The observed radio-IR luminosity ratio is 20 times 
higher than the value predicted by t he radio-infrared re- 
lation derived bv lPierini et al.l ()2003[ ). given the observed 
amount of i nfrared luminosit y. This reinforces the result 
obtained bv IXu et al.l (|2003[ ). suggesting that the accel- 
eration of relativistic particles in the shock region is not 
associated with star formation. 

4) 70% of the total star formation in SQ is associated 
with the extended FIR emission. The star formation effi- 
ciency for this component is rather higher than that ob- 
served in spiral galaxies. This could in principle be due to 
a new mode of star formation operating in the group and 
involving cooling of the X-ray halo gas, possibly enhanced 
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by dust-plasma particle collisions. However the lack of 
cold gas measurement for this distributed star formation 
component makes this hyphothesis rather uncertain. In 
addition, the evidence of tidal interactions and star for- 
mation in tidally displaced cold gas, like SQ A or SQ B, 
suggest that the distributed star formation observed in the 
group is in part occuring out of gas of galaxian origin. 

5) For both the shock region and the extended emission 
we have measured substantial 8/xm flux, mainly due to 
PAH molecule line emission and tracer of cold gas phases. 
This, together with simple estimates of the dust injection 
rate across the shock front and in the halo of the group, 
demonstrates that the mechanism powering most of the 
dust emission seen on the FIR residual maps is not dust 
coUisional heating. However, given the high FIR-X-ray 
luminosity ratio, the current data do not exclude that the 
majority of the hot gas cooling can still be due to dust- 
plasma particle collisions rather than X-ray emission. For 
dust in the group halo we estimated that dust collisional 
heating can power up to 10*^erg/s of the dust luminosity. 



6) We have found that the morphology of the FIR dust 
emission in the Seyfcrt 2 galaxy NGC 7319 does not cor- 
respond to the cold gas distribution. Therefore, the FIR 
emission in this galaxy does not seem to be related to star 
formation. Given the AGN UV luminosity, estimated from 
the observed hard X-ray flux, and from considerations on 
the particular torus-galaxy disk geometry, we deduced 
that accretion powered radiation from the AGN, propa- 
gating through the disk of the galaxy, can power most of 
the observed dust infrared emission. 
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Fig. 1.— Stephan's Quintet Multiwavelength Data. Upper row from left to right: SDDS r-band, GALEX FUV, VLA radio , XMM NEWTON soft X-ray. Lower row: SPITZER IRAC 
8/im and MIPS 24, 70, 160/im. Crosses on the maps identify the galaxy centers. The position (0,0) coincides with RA = 22'^ 36'02. 4", Dec = +33°57'46.0". The units on the axis are 
arcmin. At the distance of SQ (94 Mpc) every ar cmin corresponds to 24kpc. Note that NGC 7320 is a foreground galaxy at the distance of 10 Mpc and its HI distribution, not shown in 
this figure, can be found in [Williams et al.l ||200^ 1. 



20 



Natale et al. 




-0.32 

-1.00 
133.18 

43.84 



i 



Fig. 2. — Fit results for the lO/im map. Top-left: original lO/im map; Top-right: best fit map; Bottom-left: fit residuals; Bottom-right: 
"deconvolved" map. The crosses on the original map identify the centers of the fitted sources. Units of the values aside the color bars are 
MJy/sr. Note: The fit residual map, shown on the bottom-left panel, has all the fitted sources subtracted. It differs from the 70fim "FIR 
residual map" , whose contours arc shown in the upper panels of Fig. [J] where the SHOCK source has not been subtracted. 
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Fig. 3. — Fit results for the WOfim map. Top— left: original 160^m map; Top-right: best fit map; Bottom-left: fit residuals; Bottom-right: 
"deconvolved" map. Units of the values aside the color bars are MJy/sr. Note: The fit residual map, shown on the bottom-left panel, has all 
the fitted sources subtracted. It differs from the 160/im "FIR residual map" , whose contours are shown in the lower panels of Fig. |4] where 
the SHOCK source has not been subtracted. 
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Table 1. Results of the FIR map fitting 



Source 


RA 


Dec 




^160/:iin 




■fieO/im 




(h m s) 


( ° ' ") 


(arcs) 


(arcs) 


(mJy) 


(mJy) 


SQ A 


22 35 58.81 


33 58 51.1 


16 


15 


134 ±15 


265 ± 86 


HII SE 


22 35 59.08 


33 57 34.0 


14 


40 


61 ±9 


186 ± 140 


HII SW 


22 35 56.10 


33 57 43.1 


11 


ND 


31 ±5 


ND 


SQ B 


22 36 10.54 


33 57 20.5 


14 


36 


63 ±8 


302 ± 34 


N7319 


22 36 3.76 


33 58 33.0 


12 


30 


616 ±62 


1192 ±130 


N7319 (HII) 


22 36 3.01 


33 59 0.0 


24 


ND 


50 ±8 


ND 


N7320 


22 36 3.38 


33 56 52.7 


44 


48 


716 ± 73 


1899 ± 196 


N7320 (HII) 


22 36 4.54 


33 57 4.5 


16 


ND 


111 ± 13 


ND 


HII N 


22 36 3.13 


33 59 34.6 


23 


42 


52 ±6 


98 ±41 


SHOCK^ 


22 36 0.10 


33 58 0.3 


61 


58 


258 ± 30 


1218 ±230 



Note. — Col. 1: FIR source; cols. 2-3: Source center RA and Dec coordinates; 

cols. 4-5: Averaged FWHMs of the sources at 70/im and 160/im; cols. 6-7:Total flux 
densities at lOi^m and 160idin. Some of the sources are not detected at IdO/^m (ND = 
"non detected") 

''The SHOCK source on this tabic refers only to the gaussian used to fit the emission 
centered in the middle of the shock region. Therefore the quoted fluxes do not represent 
the actual fluxes coming from that area. The emission from the shock region is given 
in the Table 2. 



Fig. 4. — FIR residual maps (see Sect. 13. 2p countours overlaid on HI (left panels), X-ray (middle panels) and FUV maps (right panels). The countours in 
the upper row are from the 70/xm residual map while those in the lower row are from the 160/xm residual map. The 70^im residual map contour levels are 
-0.3,0.3,0.6,0.9,1.2,1.5,1.8,2.1,2.25 MJy/sr and those for the 160^m residual map are -0.4,0.4,1.2,2.4,3.6,4.8,6.0,6.6,7.2 MJy/sr (the dashed lines are the nega- 
tive contours, corresponding to —la level). The field of view is different in the three maps. The coordinates corresponding to (0,0) are for the HI map RA = 22''36'04.1" 
and Dec = -(-33°57'47.2", for the X-ray map RA = 22''35'59.7" and Dec = -h33°58'05.9", for the FUV map RA = 22''36'01.6" and Dec = -f 33°57'42.l". The axis units 
are arcminutes for the HI and FUV maps and arcseconds for the X-ray map. Note that point sources have been removed from the X-ray map (see Fig. [l|. 
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2. Inferred infrared fluxes for the discrete sources detected on the Spitzer 70/ini map 



Source 












(mJy) 


(mJy) 


(mJy) 


(mJy) 


SQ A 


7.9 ±0.5 


11 ±2 


134 ± 15 


265 ± 86 


HII SE 


3.4 ±0.4 


7.5 ± 1 


61 ±9 


186(51) ± 142(40) 


HII SW 


2.5 ±0.1 


2.4 ±0.35 


31 ±5 


ND 


SQ B 


3.1 ±0.15 


5.6 ±0.8 


63 ±8 


302(120) ±34(13) 


NGC 7319 


68 ±3 


185 ±8 


666 ± 67 


1192 ± 130 


NGC 7320 


45 ±3 


38 ±2 


827 ±85 


1899 ± 196 


HII N 


0.69 ±0.07 


1.1 ±0.2 


52 ±6 


98(78) ±41(33) 


Shock region 


8.1 ±0.8 


6± 1.5 


80 ±30 


506 ± 200 


Extended emission 


29 ±2 


40.6 ±5 


233 ± 30 


805 ± 200 



Note. — The 8/im and 24/im fluxes are obtained from aperture photometry 
while the 70/.tm and 160/im fluxes for the galaxies and compact star formation 
regions are the results of the FIR map fitting technique. The FIR fluxes from 
the shock region are obtained by the fitting of the FIR residual maps. The MIR 
and FIR values for the extended emission are derived from the curves of growth 
of the MIR and FIR maps (see Sect]?]). The values in brackets are the 160/im 
fluxes within the size of the 70/im emission. 
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Fig. 5. — "Deconvolved" 70/Lim and 160/im SQ A emission contours overlaid on 8/im (upper row) and 24/im (lower row) maps. Units on 
the axis are arcseconds. 
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Fig. 6. — Apertures for the photometry of each souree at 8/xm. Integration areas are dehmited in red while background areas in green. The 
straight hne passing through the AGN galaxy and SQB is a map artifact. 
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Fig. 7.— 



Apertures for the photometry of each source at 24/im. Integration areas are dehmited in red while background areas in green. 
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Fig. 8. — Ratio between fluxes inferred from the fitting technique and from aperture photometry for the sources on the 24^m map. 




SQ dust emission and star formation 



29 



1x10' 



1x10- 



1x10^ ^ 



1x10' 



1x10° ^ 



1x10" 



"T 1 1 1 1 1 1 1 r- 



"T 1 1 1 1 1 r- 




1 60fim 



70/i.m 
24yLtm 



FUV 



1 00 , 150 
Rodius [arcsj 



200 



250 



Fig. 10. — Radial curves of growth of the FUV, MIR and FIR emission after masking or subtraction of star formation regions and 
the two galaxies NGC 7320 and NGC 7319 (see Sect l4.3ll . The zero on the X-axis corresponds to the center of the shock region: RA = 
22^35'59.7,nec = +33°58'14.2" 
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A [fJ-m] X [fJ,m] 

Fig. 12. — Diffuse photon heated dust SEDs for different choices of Xisrf ^^<i Xcolor- In the left panel the SEDs are calculated assuming 
Xcolor = 1 a-iid several values of Xisrf (tlic higher curves correspond to higher values of Xisrf)- In the right panel we assumed Xisrf = 1 and 
different values of Xcolor (similarly, the higher curves correspond to higher values of Xcolor)- The SEDs are calculated assuming a distance 
equal to 94 Mpc (SQ distance) and a dust mass = 10^ Mq. 




Fit;. 13. — Theoretical coUisionally heated dust SEDs for several values of the size distribution exponent k. From the lowest to the highest 
curve the k values are respectively 1,1.5,2,2.5,3,3.5,4,4.5. The SEDs are calculated assuming for the plasma density and temperature 
respectively uh = 0-01 cm~^ and T = 3 x 10® K, a dust mass M<j = lO'^ Mq and a distance equal to 94 Mpc. 
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Fig. 14. — SEDs of discrete SQ star formation regions; the plotted curves are the best fit SED (solid line), the contribution from PDR/HII 
regions (dotted curve) and diffuse dust (dashed line). The values Ks^m, K24ij,in, KTOij,in, J^"l60fMn, shown in each plot of this and the following 
figures, are the color corrections applied to each fitted point. 
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Table 3. Results from the infrared SED fitting 



Source 




Xisrf Xcolor 


F24 


J abs, local 


T abs, diffuse 


J abs, diffuse 






(10^ Mq) 






(10*3 £18) 


(10*3 S£g) 


(10*3 siS) 


(10*3 HIS) 



SQ A 


6+"-^ 

U.U_g 2 


2 


2 


0.65 (HII) 


0.45 


0.61 


0.36 


1 4+0-2 


HII SE 7318b 


n 22+° "* 


2 


0.25 


0.86(HII) 


0.42 


0.10 


0.12 


64+0 00 


HII SW 7318b 


11+001 

U-ii_o.ol 


4 


1 


0.50 (HII) 


0.08 


0.14 


0.12 


34+0 03 


SQ B 


„ ,-+0.09 
^■■"-0.04 


1 


4 


0.80 (HII) 


0.31 


0.28 


0.10 


70+0 02 


NGC 7319 




2 


4 


0.90 (AGN) 


4.13 


4.01 


1.47 


9 6+°-3 


NGC 7320 




1 


4 


0.45 (HII) 


0.014 


0.059 


0.021 


094+0 005 


Shock region 


4.4+0 = 


0.3 


4 


0.35 (HII) 


0.14 


1.00 


0.37 




Extended emission 


47^3 


0.1 


1 


0.71(HII) 


1.71 


1.56 


1.30 


4 6+°-^ 



Note. — Col. 1: Fitted source; col. 2: Cold component dust mass; cols. 3-4: Radiation field Xiarf and Xcoior parameters 
(see Sect. 5); col. 5: Fraction of 24/^m flux contributed by HII regions or AGN torus templates; col. 6: Luminosity of dust 
in PDR/HII regions or in AGN torus in the case of NGG 7319; col. 7: Luminosity of diffuse dust powered by the absorption 
of optical photons of the diffuse radiation field; col. 8: Luminosity of diffuse dust powered by the absorption of UV photons 
of the diffuse radiation field; col. 9: Total dust luminosity. 




Fig. 15. - SED of NGC 7319 performed using two different combinations of templates. In the left panel the best fit SED (solid line) 
calculated using a PDR/HII regions template (dotted curve) and a diffuse dust component (dashed line). In the right panel the best fit 
calculated using an AGN torus dust emission template (dotted) and a diffuse dust component (dashed). 
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Fig. 16. — SED of the shock region. The upper panel shows the fit performed using an PDR/HII region template plus a collisional heating 
component. The hot gas parameters adopted to calculate the coUisionally heated dust SED axe: n = 0.016 cm~^ and T = 3 x 10^ K. The 
lower panel shows the fit performed using the PDR/HII region template and a uniformly photon-heated dust model. 




Fig. 17. — SED of the extended emission. The upper panel shows the fit performed using a PDR/HII region template plus a collisional 
heating component. The hot gas parameters adopted to calculate the coUisionally heated dust SED are: n = 0.001 cm~^ and T = 6 x 10® K. 
The lower panel shows the fit performed using the HII region template and a uniformly photon-heated dust model. 
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Fig. 18.— SED of NGC 7320 and HII N: the plotted curves are the best fit SED (sohd fine), the contribution from PDR/HII regions 
(dotted curve) and diffuse dust (dashed line). 



Table 4. Measured optical-UV luminosities, star formation rates and gas amount for the detected SQ dust emitting sources 



Source 


(10*°=?) 


SFRho, 
(^) 

\ yr y 


rndircct 

(10^' 


.^>dust,UV 
(1042 




SFRxjv 
(^) 

V yr y 


SFRuv/A 
(10-3^^) 

^ yr kpc^ ' 


(^) 


S(ff2) 
(^) 


S(H+) 
(^) 




SQ A 


9. ± 2 


0.7 ±0.2 


2.75 


8.09 


4.49 


0.78 


6.4 


7.9 


6.0 




0.002 


HII SE 


6.6 ± 1.5 


0.5 ± 0.1 


2.32 


5.44 


3.02 


0.58 


5.8 


< 0.40° 








HII SW 


3.9 ± 0.5 


0.2 ± 0.04 


2.11 


1.99 


1.10 


0.35 


5.5 


3.2 


< 1.7° 






SQ B 


1.3 ± 0.3 


0.2 ± 0.05 


0.31 


4.18 


2.32 


0.28 


2.7 


5.2 


4.0 




0.002 


NGC 7319 






2.01 






> 0.22 












Shock 






6.98 


5.09 


2.83 


1.05 


3.1 




9.1 


3.0 


0.007 


Extended 






22.7 


30.21 


16.78 


4.26 


1.0 






2.4 


0.01 



Note. — Col. 1: source name; col. 2: Ha luminosity; col. 3: star formation rate derived from the Ha and 24/j.m luminosity using the relation given 
by Calzetti et al.07; col. 4: observed GALEX FUV luminosity density corrected for Galactic extinction; col. 5: UV powered infrared dust luminosity; 
col. 6: absorbed UV luminosity density derived by dividing the UV powered dust luminosity Ldiist.uv by A^'{UV) — 1.8 x 10"'"'' Hz; col. 7: SFR based 
on the inferred total (obscured and unobscured) UV luminosity density, using the relation by Salim et al.07; col. 8: star formation rates in col. 7 divided 
by the projected emission area; col. 9: average neutral hydrogen mass column density; col. 10: average molecular hydrogen mass column density; col. 
11: average X-ray emitting plasma mass column density; col. 12: dust to gas ratio (note that for each source only the corresponding available gas 
masses are considered for the estimate of Z^) 

^We measured the Ha fluxes from the interference filter maps published in Xu et al.l999. The quoted flux uncertainties are the quadratic sum of 
background fluctuations and flux calibration uncertainty (« 10%). 

^The dust to gas ratio is equal to = Md/(1.4 x Mgas) where is the dust mass inferred by the SED fit, Mgas is equal to N(H) x irR^ with 
R = FW H M (rad) X distance the source size as derived from the 70jum source fit (the factor 1.4 is applied to take into account the contribution from 
noble gasses). 

*^The value of Y^{HI) for HII SE is an upper limit because HI is not detected at that position on the maps by Williams et al.02. No value is given for 
E(i/2) because the area wasn't observed by Lisenfeld et al.2002. The value of E(i?2) for HII SW is an upper limit because CO has not been detected 
at that position. 
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X [/xm] 

Fig. 19. — SED of the total SQ dust emission. The plotted points are the sum of the fluxes of all the dust emitting sources in SQ: the 
discrete sources SQ-A, HII-SE, HII-SW, SQ-B, the AGN galaxy NGC 7319, the shock region and the extended emission (the applied color 
corrections are: 0.62 (8/xm), 0.97 (24/im), 0.90 {70/jm) and 0.97 (160/im)). The plotted curves are the total emission SED (black), the discrete 
star formation region total SED (blue), the AGN galaxy NGC 7319 SED (red) and the X-ray emitting region total SED (green). All the 
SEDs representing the total emission from more sources are derived summing the fitted SEDs as described in Sect. 5. In the case of the X-ray 
emitting region SED, we combined the SED fits for the shock region and the extended emission performed using purely photon— heated dust 
emission components (see Sect. 5.3 and 5.4). 



Table 5. Galaxy optical photometry and stellar masses 



Source 


9 


r 




NGC 7319 


13.58 


12.83 


1.6 


NGC 7318a 


15.00 


14.06 


0.8 


NGC 7318b 


13.31 


12.55 


2.1 


NGC 7317 


14.58 


13.78 


0.8 



Note. — Col. 1: source name; 
cols. 2-3: SDSS g and r apparent mag- 
nitudes; col. 4: galaxy stellar mass 
(IO^Mq) 

^The stellar mass has been estimated 
using one of the relations given by Bell 
et al. 2003 to calculate the stellar mass 
to luminosity ratio from SDSS magni- 
tudes: logio{M/Lg) ^ -0.499-f 1.519X 
{9 ~ (note that M/Lg is in solar 
units). 
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log Sh [Mq pc"^] 

Fig. 20. — Star formation rate per unit area Ssfr (M0/yr/kpc^) vs total gas mass surface density Eh (Mq/pc^) for dust emission sources 
in SQ. The plotted lines are the relations found by Kennicutt et al.07 and Biegiel et al.08. The arrows pointing towards right /left for some 
of the plotted measurements indicate that the corresponding Sh value is a lower/upper limit (see Table 4). In case of HII SE the plotted gas 
surface density corresponds to the upper limit of the atomic hydrogen column density but the total gas density can also be higher since the 
amount of molecular gas is unknown. In the case of the shock region point, whose Sh value is the sum of both cold and hot gas components, 
a bar extends towards the left until the point where Sh corresponds to the cold gas mass surface density for this source. 
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APPENDIX 

MULTISOURCE FIT OF SQ FIR MAPS 

The technique we developed for the fitting of the SPITZER FIR maps is similar to the CLEAN algorithms used in 
radio astronomy but with some additional characteristics: 1) all the sources are fitted at the same time; 2) an intrinsic 
finite source width is allowed. Our basic assumption is that all the sources we model are sufficiently well described by 
elliptical gaussians or a combination of them. This is not necessarily true because some sources can, in principle, differ 
significantly from this simple shape. However this simple approach has been sufficient for a good fitting of the FIR maps. 
The function F[x,y)^ describing an elliptical gaussian on a map, is defined by the following formulas: 

F(a;, y) = Sloe + A X (Al) 

U = {x'/a,f + {y'/ayf (A2) 

x' = [x — Xg) cos — {y — Us) sin 9 (A3) 

y' = {x - Xs)iiivLe + {y ~ ys)cos9 (A4) 

The first two equations define the shape of the elliptical gaussian while the last two are the coordinate transformations 
between the map and the frame defined by the elliptical gaussian axis. For each source the fitting procedure has to provide 
7 parameters: A (the gaussian amplitude), {xs,ys) (the position of the source center), (axiC^y) (the gaussian widths in 
two ortogonal directions), (the rotation angle of the gaussian axis from the array axis) and an offset value i?ioc (the 
local background for each source). 

A real image is not just the sampling of the original signal, of course, but it's the sampling of the convolution of the sig nal 
with the point spread function (PSF) of the telescope optics-detector instrument. Using the program STINYTIIVP^ we 
obtained theoretical PSFs of the MIPS instrument for both 70 and 160 /im bands. Fig. IE21l and lE22l show the theoretical 
PSFs for the 70 /xm and 160/im bands and their profiles. In panel (b) of these figures, it is shown also the radial profile 
of an empirical PSFs (kindly provided by G. Bendo). Although there are some intrinsic differences between empirical 
and theoretical profiles, we decided to use theoretical PSFs because they can be sampled at any desired rate and do not 
contain noise that would be added to the fit (however we made some tests to determine the level of uncertainty introduced 
by using theoretical PSFs instead of real ones in our fitting procedure. We have found that the difference on the final 
results is negligible compared to other sources of error, e.g. only « 5% on the inferred source fluxes). 

The first step of our fitting procedure is to assume a distribution of elliptical gaussian sources, with initial trial values 
of the parameters, in an array of the same size of the data. First from a direct inspection of the FIR maps and then 
after some iterations, we identified 10 main sources to be fitted using elliptical gaussians (see Fig. [2]): 2 for NGC 7320 
(one for the disk emission plus one for a compact HII region), 1 each for SQ A and SQ B, 2 for NGC 7319 (one for 
the circumnuclear emission plus one for a HII region), 1 for a source located north of NGC 7319 (HII N), 2 for NGC 
7318b (the 2 compact HII regions located on the southern spiral arms HII-SE and HII-SW), 1 for a source peaked in 
the middle of the shock region. Convolving this array with a theoretical PSF, an artificial image is obtained that can be 
compared with the real data. Varying the parameters of the elliptical gaussians, it is possible to improve the agreement 
between artificial and real images until a certain degree of accuracy is reached. To mini mize and find the best fit model 
parameters, we used the Levemberg-Marquardt (LM) algorithm fsee iBevington |[l992[ ). implemented by the IDL routine 
MPFIT2DFUP^3. Although this algorithm is less dependent on the initial trial values of the parameters, compared with 
other minimization methods, the results are inevitably dependent on these. Therefore it is a good rule to assign the 
initial values as close as possible to the real ones (or better to which we believe are the real ones). The criteria we used 
to choose the initial values for the free parameters are the following: 

1) Amplitude: any number of the order of the peak amplitude of the source (this fitted parameter has been found almost 
independent from the initial value); 

2) Position: initial positions near the peaks of the sources; 

3) Gaussian Widths: derived approximately from the apparent width of the source; 

4) Rotation angle: any number. 

5) Local Background: initial value equal to zero. 

As one can notice, the choice of the initial values does not require to know in advance any parameter but the approximate 
positions of the source centers and roughly the sizes of the gaussians. All this information can be easily derived from the 
data. The fitting procedure improves the quality of the results (in the sense that we obtain a smoother fit residual map 
and a lower reduced x^) if '^6 consider only the regions close to the sources for the fit. We created a mask to exclude 
all the pixels that are too distant from the calculation, that is more than 2-3 times the apparent sizes, from the source 
centers. 

The fitting technique as described until now has been applied successfully to the TOfim FIR map. The fitting of the 
160^m map has been performed in a slightly different way because one cannot clearly see on this map all the compact 
sources that are detected on the 70/im. This effect is most probably due to both lower resolution and intrinsic color 
differences. Since our goal was to produce consistent SEDs of the main sources in SQ, we introduced new constrains for 
the fit of the 160/zm map. We assumed that the centre position of each source on the 160/im map is the same as derived 

1^ by John Krist for the Spitzer Science Center: http://ssc.spitzer.caItech.edu/dataanalysistools/tools/contributeci/generaI/stinytim7] 
written by Craig Markwardt; |http://cow. physics. wisc.edu/~craigm/idl/fitting. html for details on the minimization routine 
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from the 70 /xm map fitting. We only allowed a common shift of all the source positions in order to correct possible small 
differences in the astrometry of the two maps. Then we assumed that the orientation angles and the ratio of the elliptical 
gaussian axial widths {axTCy) for the compact sources (SQ A, HII SE, HII SW, SQ B, HII N and the HII regions in 
NGC 7319 and NGC 7320) are the same as those inferred from the 70/im map fit. Therefore we assumed that the shape 
of the emitting compact sources is similar at 70 and 160/im even if the size can be different. This is actually expected 
because colder emission can come from regions that are simply farther away from the central heating source in case of 
star formation regions. 

Figs. [2] and [3] show the results of our fitting technique for the 70/im and 160/im maps in four panels: a) the background 
subtracted Spitzer FIR map of SQ, b) the best fit image obtained with our method, c) the fit residuals and d) the 
"deconvolved" map (that is actually the map showing the elliptical gaussians whose PSF convolution gives the best fit 
map). As one can see, this method produces best fit maps with remarkable similarity to the original maps (xfo^mZ-^froe = 
0.8, Xieopm/^frcc = 0.7). The "deconvolved" map at 70/im shows, apart from the emission of the foreground galaxy NGC 
7320 and AGN galaxy NGC 7319, many discrete compact sources that are responsible for the several peaks seen on the 
real map. There is a rather extended source that peaks in the middle of the shock region, even though some excess of 
emission is still seen on the fit residual map at this position suggesting that the accuracy of the fit is not very good for this 
faint source. At 160/im the emission is dominated by the AGN galaxy, the foreground galaxy and the source peaked on 
the shock region. It is quite interesting that the gaussian used to fit the emission at the position of the shock is very well 
aligned along the main axis of the shock ridge. However the east-west width of this gaussian is too large to be related only 
to the X-ray ridge (see Sect. 14. 2p ). Some compact star formation regions are undetected (or only marginally detected) at 
160/im, as we somehow expected because no peaks are clearly seen at their position. 

The main fitted parameters and the derived fiux densities for each of the fitted sources can be found in Table [TJ 
The uncertainties on the integrated fluxes are derived from the covariance matrix provided by the LM minimization 
algorithm. For the estimate of the error on each point of the map, necessary to evaluate we considered the variance 
of background fluctuations on each map. On the final integrated fiuxes for each source a 10% error is added due to flux 
calibration uncertainties (MIPS data handbook). As one can notice from Tabled the inferred 160/im sizes for the AGN 
galaxy NGC 7319 and the compact sources HII SE, SQ B and HII N are much larger than the 70/im size. In case of 
the AGN galaxy the size difference can be physically understood if one thinks that additional warm emission come from 
the central part of the galaxy and, therefore, the 70/xm flux is more peaked towards the center (see also discussion in 
Sect. 17. 3p . The larger 160/J,m size for the compact sources arises because of the confusion with fainter sources close to the 
main central source. For these sources we derived from the deconvolved 160/im gaussians the amount of flux within one 
70/im FWHM from the peak. These reduced 160/im fluxes are quoted within brackets in Table [2] 

Comparison between results obtained by the fitting technique and aperture photometry 

Since in Sect. [S]we performed source SED fitting including MIR and FIR measurements, performed respectively using 
aperture photometry and the map fltting technique, it is important to understand how different the fluxes obtained by 
the fitting procedure are compared to those we would have obtained by aperture photometry on a high resolution map. 
To quantify this difference, we made a test using the higher resolution 24//m map where dust emission morphology looks 
remarkably similar to the 70/xm emission (at least if one considers the main sources of emission). We convolved the 24/im 
MIR map to match the resolution of the 70/im map, resampled the map to match the pixel size of the 70/im map and then 
performed the gaussian fitting technique on the convolved map. The convolution has been done using the kernel function 
created by Carl Gordon (Gordon et al. 2008). For the fit we assumed the same distribution of gaussians as in the fit of 
the FIR maps, but leaving all the parameters free. The convolved map and the fit results are shown in Fig. IE23I In Table 
IE6I and in Fig. [S]we compared the results from aperture photometry (described in Sect. 14. ip and gaussian fitting. For all 
but one source (HII SW) the fluxes derived using the two different methods are consistent within the uncertainties and 
also the inferred source sizes are close to those found at 70/im (the problem with HII SW arises because of the vicinity of 
the nucleus of NGC 7318a that is much brighter in the MIR than in the FIR). Based on the results of this test, we are 
confident that the two photometric techniques give consistent results. 

DUST SED TEMPLATES AND MODELS 

In order to fit the observed source SEDs, we needed to model the dust emission from several environments: dust in 
PDR/HII regions, dust heated by diffuse radiation fields, dust in AGN tori and collisionally heated dust embedded in hot 
plasmas. In this section, we describe the characteristics of the SEDs we used to reproduce the emission from each of these 
regions. 

Dust in PDR/HII regions 

Dust in HII and photodissociation regions (PDR), close to young stars, is warm and emits predominantly in the MIR. To 
model this emission, we used an SED template de rived by the fitting of Milky Way star formation region infrared emission 
with the t heoretical ni o del of iGroves et al.l (|2008[ ). The Galactic star formation regions considered are the radio-selected 
sample of IConti et al.l (|2004[ ). The template is the average fitted spectra of the observed emission from these regions. 
Further details will be given in a forthcoming paper (Popescu et al. 2010, submitted. Sect. 2.9). The template is shown 
in the left panel of Fig. [TTJ 
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Dust heated by diffuse radiation fields 

To fit the diffuse dust emission, we created a grid of theoretical SEDs of emission from dust heated by a uniform radiation 
field. The code we used to calculate dust emission, developed by Joerg Fischera (details will be given in Fischera et al. 
2010, in preparation), assumes a dust composition including graphite, sil icates, iron and PAH molecules. The abundances 
and size distributions are the same as in Table 2 of iFischera etahl ()2008[ ). Dust emission is calculate d taking into account 
the st ochastic temperature fluctuations of dust grains follow ing the num erical method described by iGuhathakurta et al.l 
(|l989t ). combined with the step wise analytical solution of IVoitI (|l99lh . We created the set of SEDs by varying the 
intensity and the color of the radiation field heating the dust and keeping the dust com position fixed. The s tandard shape 
of the radiation field spectra, adopted by the dust emission code, is that derived bv iMathis et al.l (|1983( ) to model the 
Galactic local ISRF. In this spectrum, the stellar contribution consists of four components: a UV emission from early 
type stars plus three blackbody curves at temperatures T2 = 7500 K, T3 = 4000 K and T4 = 3000 K used to reproduce, 
respectively, the emission from young/old disk stars and Red Giants. Each blackbody curve is multiplied by a dilution 
factor Wi that defines its intensity. The whole spectrum is multiplied by a parameter, Xisrf , whose value is unity when the 
spectrum intensity is the same as in the local radiation field. In order to consider different colors of the radiation field, we 
defined a second parameter, Xcoior, that multiplies only the blackbody curves used to model the old disk population and 
the Red Giant emission (blackbody temperatures T3 and T4). For a range of values of Xisrf (0.1,0.3,0.5,1,2,4) and Xcoior 
(0.25,0.5,1,2,4), we calculated the corresponding dust emission and created the set of models for the SED fit. The result 
of increasing the Xisrf value is that the overall intensity of the radiation field is higher and, therefore, the equilibrium dust 
temperature is higher, leading to a warmer emission in the FIR (see left panel of Fig. [T2)) . A value of Xcoior higher than 
unity increases the relative contribution of the cold stellar emission, which peaks in the optical range. These wavelengths 
are more efficiently absorbed by big grains than by PAH molecules (and small grains). As a consequence, by varying this 
parameter, one modifies the ratio between the peak intensity at the FIR, produced by big grains and the peak intensity at 
8/im, due to PAH molecules which are mainly heated by UV photons (see right panel of Fig. [T^ . A similar effect would 
have also been produced by varying the relative abundance of PAH molecules and solid grains. 

Dust emission from AGN torus 

We selected a theore tical du s ty tor us SED between those created by Feltre et al. (2010, in preparation), using the 
model presented by Fr itz et al.l (|2006( ). In this model, the characteristic of the torus, supposed to be homogeneous, are 
defined by several parameters: opening angle, external to internal radius ratio, equatorial optical depth at 9.7/im and 
two parameters, /3 and 7, that determine the radial and the angular density profile respectively according to the formula 
p{r,9) = ar-^e-^l'=°"WI, where r is the radial coordinate and 9 is the angle between a volume element and the equatorial 
plane (note that a is determined by the torus optical depth). We selected a model with a large opening angle (140°), 
with an angular gradient in the density profile (/3 = 0, 7 = 6), external to internal radius ratio eq ual to 30 and equa tiorial 
optical depth r(9.7^m) = 6. The chosen parameters are well inside the range of values found bv lFritz et al.l (|2006D from 
the fitting of a sample of Seyfert2 galaxies. Then, we extracted the output SED in the equatorial view, according to the 
Seyfert 2 classification of NGC 7319, the AGN galaxy in SQ. The extracted SED is the template we used for fitting the 
AGN torus emission and it is shown in the right panel of Fig. [TTJ 

Dust heated by collisions in X-ray plasmas 

Finally, we needed to model the emission from collisionally heated dust. For this purpose, we used a code, created by 
J. Fischera, based on the works by [Guhathakurta et al. (1989) and Voi t, (11 991) for calculating the stochastical heating 



of the dust grains and using the results of IDraine et al.l (jl979( ) and IDwea i 19871 ) to determine the heating rates due to 
dust-plasma particle collisions. Dust composition is assumed to be a mix of graphite (47%) and silicates (53%). PAH 
molecules are not included but, since they have an extremely small destruction time scale when embedded in hot plasmas 
with T > 10^ K, their abundance is expected to be negligible. In order to perform the calculations, one needs to specify 
the physical properties of the plasma, i.e. density, temperature and metallicity. The exact value of metallicity is not 
important because dust is predominantly heated by collisions with electrons, provided mainly by ionized hydrogen and 
helium. Plasma temperature and density are derived from X-ray data (see Appendix [C]) . An importan t phenomena 
that should be taken into account while preparing the SED models is dust destruction due to sputtering (jDraine et al.l 
1197 9^. Because the efficiency of sputtering depends on the grain size, we expect that the size distributio n of grains in 
hot p lasmas is different from that of dust in the cold ISM, typicall assumed to be n{a) (x a"*^ with k = 3.5 (jMathis et al.l 
II977I ). Therefore we created a set of SEDs varying 1) the plasma physical parameters (see Table lETl) and 2) the exponent 
k of the size distribution in the range [1,3.5]. A k value lower than 3.5 means that the relative number of big grains 
is higher compared to the stan dard size distributi on. Since the destruction time scale due to dust sputtering is directly 
proportional to the grain size (jDraine et al.|[T979|) . a higher relative abundance of big grains is indeed expected in hot 
plasmas. Once the heating rates due to collisions are fixed, a change in k modifies the color of the emitted radiation (see 
Fig. [131). This happens because big grains are generally colder than small grains, and do not provide warm dust emission 
due to stochastic heating. Therefore, changing their relative abundance, one can obtain colder or warmer SEDs. 

X-RAY PLASMA TEMPERATURE AND DENSITY IN SQ SHOCK AND HALO 

As we mentioned in IB. 0.51 the hot plasma physical parameters, needed for the calculation of the emission SED from 
collisionally heated dust, in both the shock and halo gas, are derived from the X-ray emission. The temperatures 
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characteristic of the plasma have been obtained bv lTrinchieri et al] (j2005[ ) fitting XMM-NETWON X-ray data and quoted 
in their Table 3. They fitted the source X-ray spectra with two temperature plasma models. The fitted colder component 
is predominant for the shock region and has a temperature T « 3 x 10^ K that we assumed to calculate the dust emission. 
For the halo region, the two plasma components have similar luminosities. In this case, we assumed the average value 
T = 6 X 10^ K for the dust emission calculation (the fitted temperatures are Ti = 3.5 x 10^ K and T2 = 8.8 x 10^ K). 
The plasma densities have been derived from the X-ray luminosity, which can be expressed as: 

L = ntn,K{T)V (CI) 

where and rie are the ion and the electron number densities respectively, A(r) is the cooling function and V is the 
volume of the X-ray emitting plasma. From a given L, A( T) and V , one can cal culate the gas density n = (ntn,,)'^'^ by 
inverting ea lCll Given t he plasma temperatures found bv lTrinchieri et al] ()2005l ). we derived A(T) by interpolating the 
CIE curve calculated bv [Sutherland et all (jl993D at two different metallicities corresponding to 100% and 30% the solar 
metallicity. Then we assumed for the shock gas an emitting volume equal to the projected shock ridge area {A — 330 kpc^) 
multiplied by a parametrized line of sight depth I, and for the halo we adopted an emitting s phere of radius 40 kpc , 
corresponding to 90" at 94 Mpc distance. The luminosity values we used are those inferred by iTrinchieri et al.l ()2005() 
j-^^shock 2.4 X 10'*^ erg/s and L^^° f« 2.3 x 10"*-^ erg/s), properly scaled to the different assumed distance (85 Mpc in that 
paper vs 94 Mpc here) . The results for different choices of the emitting volumes and metallicity are shown in Table IE71 
For the shock region plasma we obtained densities in the range 4 — 18 x 10~^ cm~'^. As expected, the minimum density is 
obtained for higher metallicity (that is, higher A) and largest volume, while the maximum density corresponds to lower 
metallicity (lower A) and smallest volume. For the halo gas, the inferred densities are 1.1 x 10^* and 1.6 x lO^"' cm^^ for 
Z = Zq and Z ~ 0.3 Zq respectively. 



IRAC AND MIPS COLOR CORRECTIONS 

From the IRAC and MIPS data handbook, the color corrections for each Spitzer band are defined as: 

/ i F,/F,o){iy/iyoy' Rdiy 
J {y/vo)^"^ Rdv 



K = ^/Sr^ {MIPS) (D2) 

/(Af)'^^^^Ai?. 

where the subscripts Vo and Xo refer to the band reference frequency or wavelength, R is the instrument spectral response 
and To = 10.000 K (these formulas derive from the conventions used to calculate the quoted fluxes on the Spitzer maps). 

ESTIMATE OF AN UPPER LIMIT TO THE DUST INJECTION RATE FROM HALO STARS IN SQ 

An upper limit to th e dust injection ra te from halo stars in SQ can be inferred from the halo R-band surface brightness, 
S'r = 24.4 mag/arcs^ ()Moles et al.lll998l) . and using theoretical predictions for the Star dust injection. We assu med that 
all the R-band halo luminosity is produced by 3 Mq stars. According to table A.l in IZhukovska et all (|2008[ ). a 3 Mq 
star with metallicity Z = 0.02 injects ~ 1.2 x 10~^ Mq of dust during its all life. Given the lifetime of such star, 
Tiifc w 3 X 10* yr, one can derived the average dust injection rate: Md = -^d/Tiife = 3.5 x 10~^^ Mo/yr. Since a star 
spends most of its life time in the main sequence phase, we approximate the R-band flux of a single halo star as that of 
a blackbody sphere with parameters characteristic of a main sequence 3 Mq star: 

Fr = vrBR(T) J (El) 

where B^{T) is the blackbody R-band flux density at temperature T = 12000 K, R = 2Rq is the star radius and 
c? = 94 Mpc is the distance from SQ. The prediction for the R-band flux density produced by a single halo star and 
observed from the earth is Fr = 1.35 x 10"^'^ ergs/cm^/s/A. Dividing Md to Fr, we obtain the dust injection rate 
per observed R-band flux density: M^/F^ = 2.6 x 10"^^ M0/yr(ergs/cm^/s/A)~^. Combining the observed R-band 
surface brightness Sk and the solid angle covering the halo region of SQ, one can obtain the total received flux 
F^^^ — 2.9 X 10"'^'' ergs/cm^/s/A. Then the dust injection rate from halo stars is equal to: 



Md^i^] F^'^' (E2) 




Using this simple approach we found — 0.075MQ/yr. This value should be considered as an upper limit since on e 
can show, using the theoretical predictions for low-intermediate mass star dust injection rates in IZhukovska et ahl (|2008[ ). 
that 3 Mq stars have the highest value of the dust injection rate per R-band flux Md/FR . 
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Fig. E22. — a) Theoretical 160/im PSF; b) Average radial profiles empirical and theoretical 160/im PSF. 
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FlK. E23. — Convolved 24^m map fit. Top-left: convolved 24/im map; Top— right: best fit map; Bottom— left: fit residuals; Bottom— right: 

deconvolved gaussian map. Units of the values aside the color bars are MJy/sr. 

Table E6. Comparison of photometry results at 24/xm obtained with aperture photometry and the gaussian fitting 

technique 



Source 


pa.p 
24um 

(mjy) 


pfit 
24um 

(mJy) 


A. fit 

(arcs) 


A tit 
(arcs) 


SQ A 


11 ± 2 


12 ± 0.5 


16 


15 


HII SE 


7.5 ± 1 


6.7 ± 0.4 


14 


12 


HII SW» 


2.4 ± 0.35 


6 ± 0.6 


11 


24 


SQ B 


5.9 ± 0.8 


5.6 ± 0.8 


15 


12 


NGC 7319'' 


185 ± 8 


195 ± 8 


12 


12 


NGC 7320*' 


38 ± 2 


40 ± 2 


44 


43 


HII N 


1.1 ± 0.2 


ND'= 


22 


ND'= 



Note. — Col. 1: Source name; col. 2: aperture photome- 
try fluxes; col. 3: fluxes from the gaussian fitting technique 
; col. 4: average source FWHMs as derived from the TOjum 
map fitting ; col. 5: average source FWHMs as derived from 
the 24/im convolved map fit. 

^This source is the only soiircc, detected by the fitting 
technique, whose inferred fiuxcs arc considerably different. 
The size of the gaussian that fit the emission on the 24/i.m 
convolved map is much larger than that found at 70/im. The 
reason for this is the contamination from the nearby nucleus 
of NGC 7318a that emits strongly in the MIR but not in the 
FIR. 

^For NGC 7319 (the AGN galaxy) and NGC 7320 (the 
foreground galaxy) the apertures are chosen such to cover the 
entire emission of the galaxies. In the fitting procedure these 
sources are modeled with two gaussians, one for the centrally 
peaked emission and one for a periphcrical HII region. The 
sizes shown in column 4 and 5 refer to the central source 
that contribute most of the flux. 



*^ND; non detected 
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Table E7. X-ray gas parameters for the shock and the halo components 



Source 


Z 


a 


los/ R n JVfgas 




(^0) 


erg cm^ ^ 


(kpc) (cm-^) (10*' Mq) 



SHOCK 


1 


3.82 


5 


0.012 


0.65 


SHOCK 


1 


3.82 


10 


0.009 


0.92 


SHOCK 


1 


3.82 


20 


0.006 


1.30 


SHOCK 


1 


3.82 


40 


0.004 


1.84 


SHOCK 


0.3 


1.92 


5 


0.018 


0.92 


SHOCK 


0.3 


1.92 


10 


0.012 


1.30 


SHOCK 


0.3 


1.92 


20 


0.009 


1.83 


SHOCK 


0.3 


1.92 


40 


0.006 


2.60 


HALO 


1 


2.54 


40 


0.0011 


10.4 


HALO 


0.3 


1.39 


40 


0.0016 


14.4 



Note. — Col. 1: Source; col. 2: assumed metallicity; col. 3: cooling 
rate; col. 4: line of sight depth (for the shock), or sphere radius (for the 
halo); col. 5: gas number density defined as n = {ntrie)^^^', col. 6: gas 
mass. 



